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The genetic code underlying protein synthesis is a canonical example 
of a degenerate biological system. Degeneracies in physical and bio- 
logical systems can be lifted by external perturbations thus allowing 
degenerate systems to exhibit a wide range of behaviors. Here we 
show that the degeneracy of the genetic code is lifted by environ- 
mental perturbations to regulate protein levels in living cells. By 
measuring protein synthesis rates from a synthetic reporter library in 
Escherichia coli, we find that environmental perturbations, such as 
reduction of cognate amino acid supply, lift the degeneracy of the 
genetic code by splitting codon families into a hierarchy of robust 
and sensitive synonymous codons. Rates of protein synthesis asso- 
ciated with robust codons are up to hundred-fold higher than those 
associated with sensitive codons under these conditions. We find 
that the observed hierarchy between synonymous codons is not de- 
termined by usual rules associated with tRNA abundance and codon 
usage. Rather, competition among tRNA isoacceptors for aminoa- 
cylation underlies the robustness of protein synthesis. Remarkably, 
the hierarchy established using the synthetic library also explains the 
measured robustness of synthesis for endogenous proteins in E. coli. 
We further found that the same hierarchy is reflected in the fit- 
ness cost of synonymous mutations in amino acid biosynthesis genes 
and in the transcriptional control of sigma factor genes. Our study 
reveals that the degeneracy of the genetic code can be lifted by envi- 
ronmental perturbations, and it suggests that organisms can exploit 
degeneracy lifting as a general strategy to adapt protein synthesis 
to their environment. 

Abbreviations: AA, amino acid; CRI, codon robustness index 

Degeneracy, the occurrence of distinct states that share 
a common function, is a ubiquitous property of physi- 
cal and biological systems [1, 2, 3]. Examples of degenerate 
systems include atomic spectra [4], condensed matter [')], the 
nervous system [2] and the genetic code [6, 7]. Degeneracy in 
physical systems is often associated with underlying symme- 
tries [i], and in biological systems with error-minimization, 
evolvability, and robustness against perturbations [s]. Degen- 
erate states that are indistinguishable under normal condi- 
tions can exhibit distinct properties under the action of ex- 
ternal perturbations [1]. This effect, called degeneracy lifting, 
allows degenerate systems to exhibit a wide range of behaviors 
depending on the environmental context [2]. The genetic code 
governing protein synthesis is a highly degenerate system since 
18 of the 20 amino acids have multiple synonymous codons 
and 10 of the 20 amino acids are aminoacylated (charged) 
onto multiple tRNA isoacceptors. Protein synthesis rates in 
living cells respond to diverse environmental perturbations, 
which raises the question of whether any of these perturba- 
tions modulates protein levels by lifting the degeneracy of the 
genetic code. Previous experiments found that both the con- 
centration of charged tRNAs as well as the occupancy of ri- 
bosomes on synonymous codons undergo significant changes 
upon nutrient limitation [9, 10, 60]. Yet whether such environ- 
mental perturbations lift the degeneracy of the genetic code by 



modulating the expression level of proteins is unknown. Here, 
we propose to use amino acid limitation in the bacterium Es- 
cherichia coll as a model system to investigate whether the 
degeneracy of the genetic code can be lifted by environmen- 
tal perturbations, and how degeneracy lifting could provide a 
general strategy to adapt protein synthesis to environmental 
changes. 



Results 

Degeneracy lifting upon amino acid limitation. We considered 
synonymous codons for seven amino acids: Leu, Arg, Ser, Pro, 
He, Gin, and Phe. This set of seven amino acids is representa- 
tive of the degeneracy of the genetic code, in that it includes 
six-, four-, three- and two- fold degenerate codon families. We 
constructed a library of 29 yellow fluorescent protein (yfp) 
gene variants, each of which had between six and eight syn- 
onymous mutations for one of the seven amino acids (Fig. lA). 
In this library, we designed each yfp variant to characterize the 
effect of one specific codon on protein synthesis. We expressed 
the yfp variants constitutively at low gene dosage (2 copies / 
chromosome. Fig. IB) in E. coli strains that were auxotrophic 
for one or more amino acids. We monitored growth and YFP 
synthesis in these strains during amino acid-rich growth as 
well as during limitation for each of the seven amino acids 
(Methods). 

During amino acid-rich growth, our measurements re- 
vealed that protein synthesis rates were highly similar across 
yfp variants, with less than 1.4-fold variation within all codon 
families (Fig. ID, grey bars). Thus, under rich conditions, the 
degeneracy of the genetic code remains intact with respect 
to protein synthesis. Strikingly, under amino acid-limited 
growth, codon families split into a hierarchy of YFP synthesis 
rates (Fig. IC, ID). We found that some synonymous codons, 
such as CTA for leucine, were highly sensitive to environmen- 
tal perturbation, causing YFP synthesis rates to be near zero 
in response to the limitation of these codons' cognate amino 
acids. Conversely, other synonymous codons, such as CTG 
for leucine, were more robust to the same perturbation with 
synthesis rates of YFP up to 100-fold higher than the sensitive 
ones^. In addition to fluorescence, this difference in robust- 
ness was reflected in protein levels measured with Western 
blotting (Fig. SI). Notably, even a single substitution to a 
perturbation-sensitive codon in the yfp coding sequence re- 
sulted in more than a 2-fold difference in YFP synthesis rate 
during limitation for the cognate amino acid, without any 
effect on synthesis rate during amino acid-rich growth (Fig. 



^We define codons as robust when the synthesis rate from the corresponding yfp variant during 
cognate amino acid limitation is higher than the average synthesis rate within that codon family. 
Similarly, we define codons as sensitive when the synthesis rate from the corresponding yfp variant 
during cognate amino acid limitation is lower than the average synthesis rate within that codon 
family. 
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S2). Only those codons that were cognate to the limiting 
amino acid caused splitting of YFP synthesis rates (Fig. S3). 
Interestingly, the splitting was more acute for codon fami- 
lies with six- fold degeneracy (Leu, Arg, Ser), while splitting 
was weaker for codon families with four-, three- and two-fold 
degeneracies (Fig. ID, first row vs. second row). These re- 
sults support the idea that greater degeneracy typically allows 
systems to exhibit wider range of responses to environmental 
perturbations [2]. In subsequent experiments, we focused on 
the two codon families, leucine and arginine that displayed 
the largest range of splitting. These two families constitute 
16% of codons across the genome of E. coU. 



Intracellular determinants of the hierarchy among synony- 
mous codons. We sought to identify the intracellular parame- 
ters that determine the observed hierarchy of degeneracy split- 
ting during amino acid limitation. To this end, we quantified 
the robustness of synthesis rate to amino acid limitation as 
the ratio of YFP synthesis rates between amino acid-limited 
and amino acid-rich growth phases. Protein synthesis rate is 
known to be correlated with codon usage and tRNA abun- 
dance during artificial over-expression of proteins [12, l:i]. 
However, we found that robustness of YFP synthesis to amino 
acid limitation was not correlated with either codon usage or 
tRNA abundance (r^ — 0.08 and 0.00, squared Spearman 
rank-correlation. Fig. S4). We then considered determinants 
of protein synthesis that might be important specifically dur- 
ing amino acid limitation. tRNA isoacceptors are uniformly 
charged (aminoacylated) at about 80% under amino acid-rich 
conditions [14, 15]. However during perturbations such as 
amino acid limitation, some tRNA isoacceptors cognate to 
the amino acid are almost fully charged while other isoaccep- 
tors in the same family have charged fractions that are close 
to zero [10, 16]. A theoretical model proposed that such se- 
lective charging arises from differences in the relative supply 
and demand for charged tRNA isoacceptors [!)]. While it is 
unclear how this mechanism could solely control protein lev- 
els, charged tRNA play an essential role as substrates for the 
elongation of ribosomes across individual codons [17]. Con- 
sequently, we hypothesized that selective charging of tRNA 
isoacceptors also underlies the observed splitting in synthesis 
rates among yfp variants. Consistent with this hypothesis, 
charged fractions of leucine and arginine tRNA isoacceptors 
during limitation of cognate amino acid starvation measured 
in a previous work [10] were correlated with the robustness of 
synthesis rates from yfp variants after accounting for codon- 
tRNA assignments (r^ = 0.78, Fig. S5). 

We experimentally tested whether varying the concentra- 
tion of charged tRNA could change the hierarchy of protein 
synthesis rates initially revealed by amino acid limitation. To 
this end, we co-expressed each one of the leucine or arginine 
tRNA isoacceptors together with each of the six leucine or 
arginine variants of yfp, respectively (Fig. 2). Previous work 
[l(i] showed that overexpression of a single tRNA isoacceptor 
cognate to a limiting amino acid enables it to compete better 
in the common charging reaction against other isoacceptors. 
As a result, charged tRNA concentration of the overexpressed 
isoacceptor increases, while charged tRNA concentrations of 
the remaining isoacceptors for that amino acid decrease or re- 
main unchanged [l(i]. We found that yfp variants constructed 
with perturbation-sensitive codons exhibited higher synthe- 
sis rates upon co-expression of tRNA isoacceptors cognate to 
those perturbation-sensitive codons (Fig. 2, A-B , bottom 
three rows, solid black-outlined squares). Conversely, yfp vari- 
ants with perturbation-robust codons exhibited lower protein 
synthesis rates upon co-expression of non-cognate tRNA isoac- 



ceptors (Fig. 2, A-B, top three rows, non-outlined squares). 
These two patterns of changes in YFP synthesis rate mir- 
ror previously measured changes in charged tRNA concentra- 
tion upon tRNA co-expression [Ki], thereby suggesting that 
the observed hierarchy in synthesis rates of yfp variants are 
tightly coupled with the concentrations of cognate charged 
tRNA isoacceptors during amino acid limitation. By con- 
trast, tRNA co-expression did not affect synthesis rates from 
yfp variants in the absence of perturbation, i.e., during amino 
acid-rich growth (Fig. 2C). We observed several codon-tRNA 
pairs with mismatches at the wobble position but that do 
not satisfy known wobble-pairing rules (Table S9), and that 
showed an increase in YFP synthesis rate upon co-expression 
of the tRNA isoacceptor during amino acid limitation (Fig. 2, 
A-B , dashed black-outlined squares). 

A codon robustness index for endogenous proteins. We in- 
vestigated whether the hierarchy of synthesis rates measured 
for the synthetic yfp variants also governs the synthesis of 
endogenous proteins of E. coli. We first devised a general pa- 
rameter, hereafter called the codon robustness index (CRI), 
to characterize the robustness of any protein's synthesis rate 
to an environmental perturbation associated with limitation 
of a specific amino acid (Fig. 3A). We defined CRI as a prod- 
uct of codon-specific weights wcodon, and we inferred these 
weights from the synthesis robustness of yfp variants to limita- 
tion for their cognate amino acid (Fig. SB). Our formulation 
of CRI is based on the simplifying assumption that each codon 
decreases protein synthesis rate by a factor wcodon that is in- 
dependent of the codon's intragenic location, the presence of 
other codons in the coding sequence, or the specific cellular 
role of the encoded protein. By definition, w codon is unity for 
codons that are not cognate to the limiting amino acid, and 
perturbation-robust codons have a higher Wcodon value than 
perturbation-sensitive codons for the limiting amino acid. 

To test the predictive power of CRI, we selected 92 E. 
coli open reading frames (ORFs) that span a broad range of 
leucine CRI values and functional categories (Fig. S7, Table 
SI). We expressed the corresponding proteins constitutively 
as N-terminus fusions with YFP^ in an E. coli strain aux- 
otrophic for leucine (Fig. 3C, Inset,). Upon leucine limita- 
tion, we found a strong correlation between the robustness of 
protein synthesis rates from the 92 ORF-jz/p fusions and their 
leucine CRI values (Fig. 3C, = 0.61, P = 10"^^, squared 
Spearman rank-correlation). Similarly, arginine CRI was also 
strongly correlated with robustness of a library of 56 ORF-y/p 
fusions during arginine limitation (r^ = 0.59, P = 10"^^ Fig. 
S8, Table S2). By contrast, standard measures of transla- 
tion efficiency under amino acid-rich conditions such as codon 
adaptation index [18], tRNA adaptation index [19] or folding 
energy of the mRNA around the start codon [20] displayed 
only a weak correlation with protein synthesis rate from the 
ORF-jz/p fusions during amino acid-rich growth (r^ — 0.10, 
0.08, and 0.02 resp.. Fig. S9). We further found that changes 
in Leu CRI calculated from the yfp data could predict both 
the effect of tRNA co-expression and that of synonymous mu- 
tations on protein synthesis from E. coli ORFs during leucine 
limitation (Fig. 3D, Fig. SIO). Importantly, similar to our 
results using yfp reporters, neither tRNA co-expression nor 
synonymous mutations for E.coli ORF-y/p fusions had a sig- 
nificant effect on the synthesis rates from these ORFs during 
leucine-rich growth in absence of environmental perturbations 



^The YFP fusion partner in the 92 ORF-y/p fusions used for testing Leu CRI was encoded by the 
CTG variant of yfp that has the highest, most robust synthesis rate during leucine limitation. Sim- 
ilarly, the YFP fusion partner in the 56 ORF-i//p fusions used for testing Arg CRI was encoded by 
the AGA variant of yfp that has the highest, most robust synthesis rate during arginine limitation. 
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(Fig. Sll). Thus the degeneracy of the genetic code under- 
hes the levels of endogenous protein production only during 
response to environmental perturbations. 

Consequences of degeneracy lifting for fitness and gene reg- 
ulation. Degeneracy splitting in physical systems can be ex- 
ploited to encode information related to the environmental 
context [21, 22]. We asked whether bacteria might similarly 
exploit the degeneracy splitting of genetic code during re- 
sponse to amino acid limitation. Hence we tested whether the 
expression of amino acid biosynthesis genes that enable bacte- 
ria to adapt to amino acid limitation is affected by the hierar- 
chy between robust and sensitive codons. We found that mu- 
tating codons that are perturbation-robust to those that are 
perturbation-sensitive in the leucine- biosynthesis genes leuA, 
leuC and leuD, and the arginine biosynthetic gene carA de- 
creased their protein synthesis rate during cognate amino acid 
limitation, but not during amino acid-rich growth (Fig. S12). 
Interestingly, in the case of leuA and carA, the same synony- 
mous mutations also resulted in a fitness cost for prototrophic 
strains upon downshift from amino acid-rich to amino acid- 
poor conditions (Fig. 4A). Thus synonymous mutations can 
have a significant fitness cost during an environmental pertur- 
bation, which is distinct from that measured under nutrient- 
rich conditions in the absence of any perturbation [20, 23]. 
However swapping codons that are perturbation-robust to 
those that are perturbation-sensitive in other biosynthesis 
genes (see argA and leuC in Fig. 4A) did not significantly 
affect fitness, suggesting that the hierarchy of robust and sen- 
sitive codons might be selectively utilized by bacteria to reg- 
ulate genes within a single metabolic pathway. 

Perturbations associated with amino acid limitation in E. 
colt can result in two distinct outcomes, depending on the en- 
vironmental conditions: On one hand, when substrates used 
in amino acid biosynthesis are still abundant in the environ- 
ment, the cell up- regulates corresponding biosynthesis genes 
to mitigate the limitation of amino acids and resume growth. 
On the other hand, in the absence of substrates for amino acid 
biosynthesis, E. coli can survive a prolonged period in amino 
acid-poor environments through a cellular response mediated 
by sigma factors [21, 25]. We found that genes encoding sev- 
eral stress-response sigma factors {rpoS, rpoE and rpoH) are 
enriched in TTA and TTG, the leucine codons that ensure ro- 
bust protein synthesis during leucine limitation (Fig. 4B, top 
panel). By contrast, genes for the housekeeping sigma factor 
(rpoD) and a few minor sigma factors {feci, fliA) are enriched 
for CTC and CTT, which are sensitive to leucine limitation. 
This contrasting pattern is observed for leucine (but not for 
arginine), and is further mirrored by the change in transcript 
abundance for sigma factor genes in response to leucine lim- 
itation (Fig. 4B, bottom panel). Hence degeneracy splitting 
in the genetic code might be exploited in concert with tran- 
scriptional control to regulate protein levels. 



Discussion 

In summary, we have found that the degeneracy of the genetic 
code does not have a role in regulating protein synthesis dur- 
ing amino acid-rich growth. By contrast, the splitting of this 
degeneracy upon reduction in amino acid supply has a potent 
effect on protein synthesis that results in up to 100-fold dif- 
ferences in protein synthesis rates between synonymous gene 
variants. Such a large role for synonymous codons in protein 
synthesis is surprising given that other post-transcriptional 
mechanisms such as protein degradation are known to play a 
significant role upon amino acid limitation [2()]. We identified 



competition between tRNA isoacceptors for aminoacylation 
as a key determinant of the hierarchy of protein synthesis 
rates during amino acid limitation. Low concentration of a 
charged tRNA isoacceptor can cause ribosomes to selectively 
pause at its cognate codon"^ and trigger ribosome jamming 
[27], translation-recoding [28], mRNA cleavage [29, 30, 31] or 
feedback-transcriptional control [•32, 33]. It will be interesting 
to find the relative contribution of these different molecular 
processes to the degeneracy lifting uncovered here*. Here, we 
have investigated the effect of a specific environmental pertur- 
bation associated with amino acid limitation in the bacterium 
E. coll. However, this type of perturbation plays a crucial 
role in the lifecycle of other bacteria such as Myxococcus xan- 
thus and Bacillus subtilis that undergo differentiation cued by 
amino acid limitation [34, 35]. Protein synthesis during such 
differentiation events might also be regulated by degeneracy 
lifting of the genetic code. Moreover, degeneracy lifting could 
be important during protein synthesis in eukaryotes, where 
clinically-important conditions such as neoplastic transforma- 
tion and drug treatment are often accompanied by a reduction 
in amino acid supply [.'5(), .'57]. Therefore lifting the degener- 
acy of the genetic code might emerge as a general strategy for 
biological systems to expand their repertoire of responses to 
environmental perturbations. 

Materials and Methods 

Summary of key methods are given below. Detailed methods for all exper- 
iments and analyses are included in the Appendix. 

Bacterial strains. All strains used in this study were obtained from 
the E.coli Genetic Stock Center (CGSC), Yale University. Different aux- 
otrophic strains were used depending on the amino acid that was limiting 
in the growth medium (Table S5). Strains were stored as 20% glycerol 
stocks at -80°C either in 1ml cryo-vials or in 96-well plates (3799, Costar). 
For experiments involving over 25 strains, a temporary 20% glycerol stock 
was stored at -20°C in 96-well PGR plates. 

Plasmids. The pZ series of plasmids [ ] were used for expression of 
all genes constructed for this study. General features of the plasmid back- 
bones are described here. Details on individual gene constructs that were 
inserted into these plasmid backbones, including DNA sequences and plas- 
mid maps, are in Appendix. A low-copy plasmid, pZS*ll \SC101* or\ (3-4 
copies/cell), AmpR (jb/a gene) and a constitutive Pj^tetOl promoter] was 
used for expression of all fluorescent reporter genes and their fusions. The 
synthetic ribosome binding site (RBS) in the original pZS*ll backbone was 
replaced by a modified T7-based RBS that resulted in efficient expression of 
most coding sequences. A medium-copy plasmid, pZA32 [pl5A ori (10-12 
copies/cell), ChIR (cat gene) and Pi^lacOl promoter] was used for expres- 
sion of all tRNA genes. Strains with pZA32 plasmids were grown with 
ImM IPTG to ensure constitutive expression of all tRNA genes. Standard 
plasmids pUClS and pUC19 were used as intermediate cloning vectors for 
site-directed mutagenesis. 

Gene synthesis and cloning. A single yfp sequence was built de novo 
(synthesis by Genscript, USA). All subsequent yfp variants were constructed 
using a site-directed mutagenesis kit (Stratagene). tRNA genes and E. coli 
ORFs were amplified from the chromosome of a wild-type E. coli strain 
(MG1655) by PGR (Details on cloning and genes sequences in Appendix). 

Amino acid limitation experiments. Overnight cultures were inocu- 
lated from glycerol stocks or fresh colonies and grown in a MOPS-based 
rich-defined medium with SOOfiAI of 19 amino acids and lOmM serine 
at 30°G with shaking. For experiments involving amino acid limitation, 
overnight cultures were diluted 1:1000 into a similar rich-defined medium 



' A recent genome-wide study found increased ribosome pausing at serine codons during serine- 
limited growth of E. coll. Interestingly, ribosomes paused significantly only at four out of the six 
serine codons, and these four codons are precisely the same ones that caused YFP synthesis rate to 
be sensitive to serine limitation in our experiments (Fig. S13). 

'^We measured the change in mRNA levels of different yfp variants in response to amino acid 
limitation. Changes in mRNA levels were correlated with corresponding changes in YFP synthesis 
rates upon amino acid limitation (Fig. S14), but were smaller than expected. 
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as the overnight cultures. However the amino acid, whose limitation was to 
be induced, was added at a reduced concentration and supplemented with 
its methyl-ester analog (see Table S6 for exact concentrations). Amino 
acid methyl-esters are inefficiently metabolized analogs of the correspond- 
ing amino acids and have been previously used for steady growth of E. coli 
under amino acid limiting conditions [51, 52] (see Figs. S15 and S16 for 
the effect of methyl-ester on growth and measured robustness during amino 
acid limited growth). Slight variations in the initial concentration of either 
the limiting amino acid or its methyl-ester only results in shifting of the 
transition to a higher or lower cell density without appreciable changes in 
growth rate (see Notes SI and S2). Growth and fluorescence were quanti- 
fied using a standard 96-well plate reader integrated with a robotic system. 
Further details on growth protocols are given in Appendix. 

Analysis of cell density and fluorescence time series. Matlab R2009 
(Mathworks) was used for all analyses unless otherwise mentioned. All 
correlations and P-values reported in this work were calculated using the 
Matlab command 'corr' with the 'Type' option set to either 'Spearman' 
or 'Pearson' as appropriate. Growth and fluorescence time series were fit 
with exponential and linear curves in the amino acid rich and amino acid 
limited growth regimes, respectively, and the onset time of amino acid 
limited growth was automatically inferred from their intersection. Protein 
synthesis rate, S was calculated as: 



Protein synthesis rate S - 



d(Fluorescence) 



Absorbance 



ci(time) 



Fig. 3B were calculated using the robustness of protein synthesis of the 
corresponding yfp variants during cognate amino acid limitation (Fig. ID). 
Based on our non-cognate limitation experiment (Fig. S2), the uui values 
for all codons other than those cognate to the limiting amino acid are set 
to be equal to 1. For illustration, we demonstrate the calculation of Wi for 
the six leucine codons below. The exact same procedure was followed for 
other synonymous codon families. Taking log2 Wi = Wi for each codon, 
and log2 (robustness during amino acid limited growth) = SR for each 
yfp variant. 



7 X WcTA 


-1- 15 X WcTG 


= SRyfp^CTA 


[2] 


7 X WcTC 


-1- 15 X WcTG 


= SRyfp^CTC 


[3] 




22 X WcTG 


= SRyfp^CTG 


[4] 


7 X WcTT 


-1- 15 X WcTG 


= SRyfp^CTT 


[5] 


7 X Wtta 


+ 15 X WcTG 


= SRyfp^TTA 


[6] 


7 X Wttg 


+ 15 X WcTG 


= SRyfp^TTG 


[7] 



[1] 



The multiplicative factors on the LHS in front of Wi correspond to the 
frequency of the Leu codon i in the corresponding Leu variant of yfp (see 
Fig. lA). The RHS is the measured (log2) robustness of protein synthesis 
from the corresponding yfp variant during Leu limitation (see Fig. ID). 
These equations were solved simultaneously to determine the Wi value for 
each Leu codon. Revised Wi values based on yfp measurements in the 
presence of *^^*^Leu2 tRNA (Fig. 2) were used for calculation of Leu CRI 
in the case of GAGLeu2 tRNA co-expression with E. coli ORFs (Fig. 3D). 



First, the above formula was evaluated at the onset time of amino acid 
limited growth using the exponential fits for absorbance and fluorescence 
data in the amino acid rich growth regime. Next, the same formula was 
evaluated at the onset time using the linear fits for absorbance and fluo- 
rescence data in the amino acid limited growth regime. These two values 
correspond to the protein synthesis rates reported for the amino acid rich 
and amino acid limited growth regimes (such as the data in Fig. ID). 
Further details of this analysis are given in Appendix. 

Calculation of CRI. CRI for a protein coding sequence corresponding to a 
limiting amino acid was calculated by multiplying the wi values for codons 
cognate to the limiting amino acid in that sequence. Wi values shown in 
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Fig. 1: Degeneracy lifting associated with amino acid limitation. 

(A) A library of 29 variants of the yellow fluorescent protein gene {yfp) was synthesized. In this library, each variant (represented 
as a horizontal line) was designed to measure the effect of one specific codon on protein synthesis rate. The identity of this 
codon and that of its cognate amino acid is indicated to the left of each yfp variant, and the locations of this codon along yfp 
are represented as thick vertical bars. Other codons for the same amino acid that were identical across all yfp variants in each 
codon family are represented as thin vertical bars. 

(B) Each yfp variant was constitutively expressed from a low-copy vector [SClOl* ori, 2 copies / chromosome) in E. coli 
strains that were auxotrophic for one or more of seven amino acids. 

(C) To induce amino acid limited growth, we adjusted the initial concentration of an amino acid in the growth medium to 
a level below that is required for reaching saturating cell density. A methyl-ester analog of the amino acid supported steady 
growth in the amino-acid limited phase. Growth and fluorescence curves for two yfp variants, CTA, red, and CTG, black, are 
shown as illustrative examples of degeneracy splitting upon limitation for the cognate amino acid, leucine. 

(D) YFP synthesis rates during limitation for cognate amino acid - blue; YFP synthesis rates during amino acid-rich growth 
- grey. YFP synthesis rate was defined as the rate of fluorescence change divided by the cell density. Synthesis rates were 
normalized by the maximum value within each synonymous codon family, and separately in the amino acid-rich and amino 
acid-limited growth phases. Normalization factors (amino acid - rich, limited): Leu - 94, 81; Arg - 89, 113; Ser - 217, 343; 
Pro - 306, 49; He - 295, 45; Gin - 185, 83; Phe - 311, 20; (arbitrary units). Error bars show standard error over three replicate 
cultures. 
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Fig. 2: Altering the hierarchy of degeneracy splitting among synonymous codons. 

The five leucine (arginine) tRNA isoacceptors were co-expressed together with each of the six leucine (arginine) yfp variants 
resulting in thirty iRNA-yfp combinations for leucine (arginine). 

(A, B) Each square in the left (right) table corresponds to the dilference in YFP synthesis rates of each yfp variant between 
the tRNA co-expressed strain and the parent strain without extra tRNA during leucine (arginine) limitation. YFP synthesis 
rates were defined in the same manner and normalized by the same factor as in Fig. ID. YFP synthesis rate of the parent strain 
without extra tRNA during amino acid limitation is shown on the left of each table (same data as in Fig. ID). tRNA isoacceptor 
names are preceded by their unmodified anticodon sequences. Solid black-outlined squares correspond to codon-tRNA pairs 
that satisfy wobble-pairing rules after accounting for known post-transcriptional tRNA modifications (Table S9). Dashed black- 
outlined squares correspond to codon-tRNA pairs that do not satisfy known wobble-pairing rules but that show a significant 
increase in YFP synthesis rate upon co-expression of the tRNA isoacceptor. ^'^^Arg2m is a non-native arginine tRNA that 
was created by mutating the anticodon sequence of the ^*^'^Arg2 gene. Standard error was less than 0.05 for all squares. 
(C) Histogram of differences in YFP synthesis rate of yfp variants upon tRNA co-expression. Amino acid limited growth: 42% 
median difference; Amino acid-rich growth: 9% median difference (n = 60, aggregated for leucine and arginine). Change in 
YFP synthesis rate between each tRNA co-expressed strain and its parent strain expressing no extra tRNA was calculated as 
a percentage of the largest value between the two YFP synthesis rates. 
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Fig. 3: Degeneracy lifting for endogenous proteins. 

(A) The effect of each codon on the synthesis rate, S, of a protein during amino acid Umitation was modeled by a codon-specific 
weight, Wcodon- The codon robustness index (CRI) for any protein coding sequence was defined as the product of wcodon values 
for all codons in that sequence that are cognate to the limiting amino acid. 

(B) Wcodon values for leucine and arginine codons during limitation for their cognate amino acids were estimated from protein 
synthesis rates of the corresponding yfp variants (Methods). Wcodon values for all codons not cognate to the limiting acid were 
set to 1. 

(C) Ninety-two open reading frames (ORFs) from the E. coli genome were cloned as N-terminal fusions to YFP downstream 
a constitutive promoter into a low-copy vector (Inset, Methods). Robustness to leucine limitation is quantified as the ratio 
of protein synthesis rates between leucine-limited and leucine-rich growth phases. This measured robustness was correlated 
with estimated Leu CRI values for the 92 ORF-y/p fusions (r^ = 0.61, squared Spearman ran k-correlation, P = 10"^''). 11 
ORFs had measured robustness below the lower limit of the vertical axis (Table SI), but were included in the calculation of 

. Protein synthesis rates were normalized by the synthesis rate for the CTG variant of yfp. Error bars show standard error 
over three replicate cultures. 

(D) Two sets of ORF-y/p fusions (21 total ORFs) were co-expressed with '^^'^Leu2 tRNA. Based on the yfp data (Fig. 2A), we 
estimated a higher CRI for the first set (11 ORFs) and a lower CRI for the second set (10 ORFs) upon '~'^''Leu2 co-expression 
{Left panel, Methods). Hence we predicted that the first set should show an increase in robustness of protein synthesis during 
leucine limitation while the second set should show a decrease. These predictions agreed with measured changes for 20 of the 
21 ORFs [Right panel, = 0.57, P = 10^*). Error bars show standard error over three replicate cultures. Several error bars 
are smaller than data markers. 
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Fig. 4: Fitness cost and transcriptional control reflect degeneracy lifting. 

(A) Four different prototrophic E. coli strains were created. Each of these strains had one of tlie four amino acid biosynthesis 
genes argA (Arg), car A (Arg), leuA (Leu) and leuC (Leu) replaced at the native locus by a corresponding synonymous mutant 
ORF. These mutants were designed such that three to five perturbation-robust codons in wild-type ORF were replaced by 
perturbation-sensitive codons in the mutant ORF (see Fig. S12, B). The strains were grown in medium supplemented with 
all 20 amino acids at 800/iM, and then diluted into a medium lacking either leucine (left panel) or arginine [right panel). 
Growth lag was calculated as the time taken by each strain to reach ODsoo of 0.3 relative to a reference culture of the same 
strain grown in of all 20 amino acids. Difference in growth lag between the leuA mutant and the two controls during 
leucine downshift [left panel) was 9.2 ± 2.8 min, P — . Difference in growth lag between the car A mutant and the two 
controls during arginine downshift ( right panel) was 7.8 ± f.2 min, P = 10~^. Standard errors were calculated over six 
biological replicates for each mutant. P- values were calculated using two-tailed t-test between the leuA or car A mutant and 
the corresponding controls. 

(B) [Top panel) Genes encoding sigma factors and leucine biosynthesis genes in E. coli are biased in their Leu CRl values, 
as quantified using a 2;-score that measures the normalized deviation from the expected CRI value based on genome-wide 
codon frequencies (Appendix). The most frequent leucine codon CTG was excluded in this analysis since its frequency varies 
significantly with expression level under nutrient-rich conditions [11]. (Bottom panel) Fold-change in mRNA abundance in 
response to leucine limitation for sigma factor genes and leucine biosynthesis operons was measured using RT-qPCR. Fold- 
change of the gapA gene was used for internal normalization. Error bars show standard error over triplicate qPCR measurements. 
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Appendix 
Bacterial strains 

All strains used in this study were obtained from the E.coli 
Genetic Stock Center (CGSC), Yale University. For amino 
acid limitation experiments, standard auxotrophic strains 
(Table S5) were used depending on the amino acid that was 
limiting in the growth medium, unless mentioned otherwise. 
Strain CP78 was used for experiments involving leucine and 
arginine limitation. This strain has been used extensively in 
previous amino acid limitation studies [42, 4.3, 44, 4.5, 15, 10] 
and its multiple auxotrophy makes it a convenient choice for 
experiments involving limitation for several amino acids. The 
auxotrophic strains corresponding to the remaining amino 
acids are from the Keio-knockout collection [46], and are 
the commonly used auxotrophic strains for that amino acid 
(http: //cgsc .biology .yale . edu/Auxotrophs .php). 

For the growth lag measurements in Fig. 4A, the pro- 
totrophic strain MG1655 (Table S5) was used as the wild-type 
background. This background strain was tagged with yfp or 
rfp at the attBX locus (this tagging was a remnant from ear- 
lier experiments not related to this work, and has no relevance 
to any results presented here). Site-directed mutagenesis was 
used to create the synonymous mutant coding sequences for 
leuA, leuC, leuD, car A, argA and argF using the protocol de- 
scribed in the section on gene synthesis and cloning. Then 
to insert these mutant ORFs into their native locus without 
any additional markers, a two-step strategy based on A Red- 
mediated homologous recombination [17] was used: In the 
first step, the respective wild-type ORF was replaced by a kan 
resistance gene, and in the second step the kan gene was re- 
placed by the mutant ORF without any additional markers by 
selecting on M9-glucose plates for prototrophy of the respec- 
tive amino acid. Plasmid pSIMS [48] was used as the helper 
plasmid and a previously published recombineering protocol 
[48] was used without any modifications. 

For RT-qPCR (Fig. 4B), a leucine auxotroph of MG1655 
was created by deleting the leuB gene using the A Red- 
mediated homologous recombination protocol outlined above. 
For Western blots (Fig. SI), the auxotrophic strains in Ta- 
ble S5 were further modified by insertion of the tet repressor 
gene at the attBX site using a previous method based on A 
integrase-mediated site-specific recombination [49]. The pres- 
ence of Tet repressor enabled inducible control of YFP expres- 
sion. The Western blots for leucine and arginine yfp variants 
were performed in an MG1655 auxotroph strain background 
instead of the CP78 strain. The CP78 strain has lower trans- 
formation efficiency which prevented integration of the tet re- 
pressor gene into the chromosome. Strains were stored as 20% 
glycerol stocks at -80° C either in 1ml cryo- vials or in 96- well 
plates (3799, Costar). In addition, for experiments involving 
over 25 strains, a temporary 20% glycerol stock was stored at 
-20° C in 96-weU PGR plates. 



Plasmids 

The pZ series of plasmids [19] were used for extra- 
chromosomal expression of genes. General features of the plas- 
mid backbones are described here. Specific gene constructs 
that were cloned into these backbones is described in the 
section on gene synthesis and cloning. A low-copy plasmid, 
pZS*ll [SC101*'on (3-4 copies/cell), AmpR [hla gene) and 
a constitutive PLtetOl promoter] was used for expression of 
all fiuorescent reporter genes and their fusions. The synthetic 
ribosome binding site (RBS) in the original pZS*ll backbone 
was replaced by a modified r7-based RBS that resulted in 



efficient protein expression from most coding sequences. A 
medium-copy plasmid, pZA32 [pl5A ori (10-12 copies/cell), 
GhlR {cat gene) and PLlacOl promoter] was used for expres- 
sion of all tRNA genes. Strains with pZA32 plasmids were 
grown with ImM IPTG to ensure constitutive expression of 
all tRNA genes. Standard plasmids pUC18 and pUC19 (In- 
vitrogen) were used as intermediate cloning vectors for site- 
directed mutagenesis. Plasmid pSIM5 (13) was used as the 
helper plasmid expressing the A-Red system for all chromo- 
somal modifications in this project (except for Tet repressor 
insertion mentioned in the previous section). 



Growth and fluorescence measurements 

Overnight cultures were inoculated either from freshly grown 
single colonies or, in experiments involving more than 25 
strains, from temporary glycerol stocks stored at -20C. 
Overnight cultures were grown in a modified MOPS rich- 
defined medium [oD] made with the following recipe: lOX 
MOPS rich buffer, lOX ACGU nucleobase stock and lOOX 
0.132M K2HP04 (Teknova, Cat. No. M2105) were used 
at IX final concentration as in the original recipe. In ad- 
dition, the overnight growth medium contained 0.5% glucose 
as carbon source, 10""% thiamine and 800/iM of 19 ammo 
acids and lOmM of serine. pH was adjusted to 7.4 using IM 
NaOH and appropriate selective antibiotics (50/^g/ml ampi- 
cillin and/or 20/^g/ml chloramphenicol) were added. Amino 
acids, glucose, thiamine and antibiotics were purchased from 
Sigma. 1ml overnight cultures were grown in 2ml deep 96- well 
plates (40002-014, VWR) at 30° C with shaking at 1350rpm 
(Titramax 100 shaker) for 12 to 16 hours. 

For amino acid limitation experiments, overnight cultures 
were diluted 1:1000 into 1ml of the same MOPS rich-defined 
medium as the overnight cultures. However the amino acid 
whose limitation was to be induced was added at a reduced 
concentration and supplemented with its methyl ester ana- 
log (Table S6). Amino acid methyl esters are analogs of the 
corresponding amino acids and have been previously used for 
steady growth of E. coli under amino acid limiting conditions 
[51, 52] (see Figs. S15 and S16 for the effect of methyl ester 
on growth and robustness of YFP synthesis) . Addition of the 
methyl esters results in a steady but limiting supply of the 
amino acid due to slow hydrolysis of the ester (see Note 1). 
Concentrations of the amino acid and its methyl ester were 
chosen such that the cultures consumed the limiting amino 
acid and entered amino acid-limited growth at an ODeoo of 
0.6-0.7 (corresponding to an ODeoo value of 0.2-0.25 in our 
96- well plate reader). Slight variations in the initial concen- 
tration of either the limiting amino acid or its methyl ester 
shift the transition to a higher or lower cell density without 
appreciable changes in growth rate (see Note 2). Except for 
a single limiting amino acid, the remaining 19 amino acids 
were present at the overnight culture concentrations during 
the amino acid limitation experiments. For proline limita- 
tion, no proline was necessary in the growth medium since 
proline methyl ester supported growth at the same rate as 
proline until the ODgoo reached around 0.6. 

Diluted overnight cultures were grown in 2ml deep 96- well 
plates for 3 hours at 30° C with shaking at 1350rpm (Titra- 
max 100 shaker). After this time interval, 3 aliquots of 150pil 
from each culture was pipetted into 3 wells of 3 different 
96- well plates (3799, Costar). Wallac Victor2 plate reader 
(PerkinElmer) was used to monitor cell density (absorbance 
at 600nm) and YFP synthesis (fluorescence, excitation 504nm 
and emission 540nm) . Each plate was read every 15 min using 
a robotic system (Caliper Life Sciences) and shaken in between 
readings (Variomag Teleshake shaker) for a total period of 6- 
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10 hours. Temperature of 30° C and relative humidity of 60% 
was maintained throughout the experiment. 

In the case of experiments without methyl ester (Figs. 
S15 and S16), the same protocol mentioned above was fol- 
lowed but the methyl esters were not added to the growth 
medium. For the RT-qPCR measurements shown in Fig. 4B, 
overnights cultures were diluted 1:1000 into the same medium. 
Then when the ODgoo reached 0.5, the cells were spun down 
at 3000g for 5 min and then re-suspended in the same medium 
but either with or without leucine. Total RNA was extracted 
(see protocol below) after 30 min of shaking at 30° C, 200rpm. 

For the growth lag measurements shown in Fig. 4A, 
overnight cultures of prototrophic strains were diluted 1:200 
into medium either with or without one of leucine and argi- 
nine. Growth lag was measured as the difference in time taken 
to reach ODaoo of 0.3 between two cultures of the same strain 
- one growing in the presence of either leucine or arginine and 
another growing in its absence. 



Gene synthesis and cloning 

All gene sequences constructed for this study are provided 
in the gene_sequences.fasta file. Plasmid backbone sequences 
are provided in the plasmid_sequences.genbank file. Primer 
sequences used for cloning will be provided upon request. For 
all primers, 18 to 22bp homologies without any special primer 
design criteria were sufficient for successful PGR amplification 
with Phusion High-Fidelity DNA polymerase (NEB). 

Initial yfp construct. All yfp variants used in this study were 
modified starting from a single yellow fiuorescent protein gene 
sequence (called yfpO in the sequence file and plasmid map). 
This yfpO sequence encoded the fast-maturing 'Venus' vari- 
ant of YFP [ j-'i]. All 238 codons of yfpO were chosen such 
that they were decoded by abundant tRNA isoacceptors for 
each amino acid. Such a choice of codons ensured that the 
native level of demand for each tRNA isoacceptor inside the 
cell was minimally perturbed by the low-copy expression of 
fluorescent reporter genes. The yfpO sequence was built de 
novo (synthesis by Genscript, USA). The synthesized yfpO se- 
quence was cloned between the Kpnl and Hindlll restriction 
sites of the pZS*ll plasmid vector using standard molecular- 
biology techniques (19). The plasmid map of the resulting 
construct, pZS*ll-yfpO is shown in Fig. S17. 

Synonymous variants of yfp- A subset of codons in yfpO cor- 
responding to 7 amino acids (Leu, Arg, Ser, Pro, He, Gin, Phe) 
were mutated to create the initial 29 synonymous variants of 
yfp (yfP^ ~ yfp29 in the gene_sequences.fasta file, sequences 
in the same order as shown in Fig. lA). The 4 yfp variants 
corresponding to Pro {yfpl9-yfp22) had all the Pro codons 
mutated to the most frequent GCG codon since the original 
yfpO sequence had a few CGA and GCT codons that are more 
sensitive to Pro limitation. Similarly, all the Phe codons in 
yfpO were mutated to the most abundant Phe codon TTT 
for the two Phe variants of yfpO {yfp28-yfp29). Both these 
groups of variants (6 total) had higher overall fluorescence dur- 
ing amino-acid rich conditions than the rest of the 23 variants. 
This higher fluorescence is likely due to changes in secondary 
structure near the ribosome binding region on the mRNA as 
a consequence of mutations near ATG. However, this change 
is common across all variants within the Pro and Phe syn- 
onymous codon groups and hence is not responsible for the 
differential response to cognate amino acid limitation mea- 
sured within these synonymous codon groups. 



For constructing the 29 yfp variants, yfpO from pZS*ll- 
yfpO was flrst cloned into a plJCW cloning vector between the 
Kpnl and Hindlll restriction sites. A commercial site-directed 
mutagenesis kit (Quickchange Lightening Multi, Applied 
Biosystems) was used to introduce the mutations correspond- 
ing to each of the 29 variants and the manufacturer's protocol 
was followed. The resulting variants were verified by Sanger 
sequencing and then cloned into the pZS*ll expression vector 
backbone between the Kpnl and Hindlll sites. The 22 single 
GTA variants of yfp (Fig. S2) were constructed using the same 
procedure as above. The 29 yfp variants for Western blotting 
(Fig. SI) were created using the same procedure as above, but 
with the addition of a 22 codon sequence at the 5' end that 
encoded a 3X-FLAG peptide recognized by a commercially 
available, anti-FLAG, antibody (Sigma). The 22-codon se- 
quence is: GACTACAAAGAGCATGACGGTGATTATAAA- 
GATGATGACATCGAGTAGAAGGATGACGATGAGAAG. 

tRNA expression vectors. The 5 distinct Leu tRNA isoaccep- 
tors encoded by the genes leuQ, leuU, leuW, leuX and leuZ, 
and the 4 distinct Arg tRNA isoacceptors encoded by the 
genes argV, argX, argil and argW were cloned between the 
EcoRI and HindHI sites of the pZA32 expression vector (Fig. 
S18). These genes were amplified by PGR from the chromo- 
some of E. coll strain MG1655. In addition to these native 
tRNA genes, a synthetic tRNA gene arg2m cognate to the 
GGA Arg codon was also created. Normally, the ^'-^''Arg2 
tRNA with ICG anti-codon reads the GGA codon inefficiently 
through a purine-purine wobble pairing. Expressing a syn- 
thetic tRNA with an anticodon UGG restores efficient read- 
ing of this codon and is equivalent to increasing the supply 
of the corresponding cognate tRNA isoacceptor. This syn- 
thetic tRNA isoacceptor was created from the pZA32-argV 
expression vector using overlap PGR to introduce the neces- 
sary single bp mutation in the anticodon of argV. The pZA32 
vectors with the tRNA genes were electroporated into strains 
already containing the YFP expression vectors. 

Library of E. CoH ORF-yfp fusions. 92 E. coli Open Read- 
ing Frames (ORFs) were selected for experimental validation 
of the Leu Godon Robustness Index (Leu GRI). These ORFs 
were chosen to span a wide range of predicted Leu GRI values 
and functional categories (Fig. S7 and Table SI). First, a 
modified pZS*ll-yfpO vector backbone was created in which 
the start codon of yfpO was replaced by a GGSGGS hexa- 
peptide linker sequence: GGTGGATGGGGGGGTTCT con- 
taining a BamHI restriction site. Next, the 92 ORFs (without 
the stop codon) were amplified by PGR from the chromo- 
some of E. coli strain MG1655 with 5'-KpnI and 3'-BamHI 
restriction site overhangs. These PGR fragments were cloned 
into the modified pZS*ll-yfpO vector backbone containing the 
BamHI restriction site. 13 of the 92 ORFs had either an in- 
ternal Kpnl or an internal BamHI site. In these cases, a larger 
fragment that included adjoining sections of the pZS*ll-yfpO 
vector was constructed by overlap PGR and then cloned us- 
ing other restriction sites (EcoRI or HindHI). Thus the flnal 
constructs had one of the 92 E. coU ORFs connected through 
a hexapeptide linker with yfpO. All the cloned sequences were 
verifled by PGR for inserts of right length and around 40 ORF 
constructs were verified by Sanger sequencing. Two biologi- 
cal replicates of each ORF construct were compared for their 
synthesis robustness values as measured during the amino acid 
limitation assay and these values showed a high degree of cor- 
relation (Pearson p = 0.93, Fig. S19). 

For validating the Arg codon robustness index (Arg GRI), 
56 E. coli ORFs that included a subset of the above 92 ORFs 
were chosen (Table S2). The cloning procedure was exactly 
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analogous to the above 92 ORFs but with one difTerence: the 
yfpO part of the fusion construct was replaced by a synony- 
mous variant of yfpO {yfp7) that had the Arg codon AGA in- 
stead of the CGT and CGC codons in the yfpO sequence. The 
codon AGA has the highest Wi value among the Arg codons 
(see Fig. 3B) and hence has a minimal effect on the measured 
robustness of the ORF fusions during Arg limitation. 

Co-expression of *^^^Leu2-tRNA with E. CoH ORF-yfp 
fusions. Out of the 92 E. coli OKF-yfp fusions, 21 were cho- 
sen for co-expression with the ''^'~'Leu2 tRNA that is cog- 
nate to the codons CTC and CTT. The 21 ORFs were chosen 
such that 11 of them had a lower Leu CRT prediction than 
their wild-type counterparts while the other 10 ORFs had a 
higher Leu CRI prediction than their wild-type counterparts 
(Table S3). This choice also corresponded respectively to ei- 
ther high frequency of the non-cognate TTA and TTG codons 
for Leu2 or high frequency of the cognate codons CTC 
and CTT. The strains containing the 21 ORF fusions were 
each made electro-competent and then transformed with the 
pZA32-leuU plasmid that expresses *^*'^Leu2. 

Synonymous variants of E. Coli ORF-yfp fusions. Out of 

the 92 E. coli ORF-y/p fusions, 13 were selected for creating 
synonymous mutants (Table S2). These 13 ORFs had a high 
frequency of one or both of the Leu codons, TTA or TTG and 
these codons were mutated to the Leu codon, CTC. All these 
3 codons, TTA, TTG and CTC occur at similar frequencies 
on average across the genome of E. coli.The 13 ORF-y/p fu- 
sions were amplified by PGR from the pZS*ll vectors between 
the EcoRI and Xbal restriction sites (see Fig. S17). These 
fragments were cloned between EcoRI and Xbal sites of the 
pUClQ cloning vector. A commercial site-directed mutagen- 
esis kit (Quickchange Lightening Multi, Applied Biosystems) 
was used to introduce the TTA, TTG CTC mutations. 
A unique primer was designed for each of the TTG or TTA 
codons in the 13 ORFs, and these primers encoded the CTC 
mutation. All the primers corresponding to each ORF were 
mixed and then used in the mutagenesis reaction. This pro- 
cedure resulted in mutant coding sequences with TTA, TTG 
— >■ CTC mutations at random locations. 10 colonies for each 
ORF were sequenced and each unique mutant sequence was 
then cloned into the pZS*ll expression vector. At the end of 
the procedure, a total of 63 constructs were created that each 
had between one and seven TTA, TTG — )■ CTC mutations 
(see gene_sequences.fasta file for exact sequences). 



Total RNA extraction 

Total RNA was extracted for two different experiments (Figs. 
4B, S14). Phenol-chloroform extraction method was used to 
obtain total RNA. Briefly, 3ml of cells were quickly mixed 
with 5ml of ice-cold water and harvested by centrifugation at 
3000g for lOmin. Cell pellets were re-suspended in 500/^1 of 
0.3M sodium acetate-lOmM EDTA, pH 4.8 buffer. The re- 
suspended cells were mixed with 500/il of acetate-saturated 
phenol-chloroform at pH 4.8, SO^il of 20% SDS and 500^1 of 
acid- washed glass beads (G1277, Sigma). The mixture was 
shaken in a vortexer for 20 min at 4C. The aqueous layer 
was extracted twice with acetate-saturated phenol-chloroform 
at pH 4.8 and once with chloroform. Total RNA was pre- 
cipitated with an equal volume of isopropanol and washed 
with 70% ethanol-50mM sodium acetate pH 4.8 and finally 
re-suspended in 200/il of RNase-free water. 20/il of the to- 
tal RNA was treated with DNase (EN0521, Fermentas) to 
remove residual DNA contamination (manufacturer's instruc- 



tions were followed). The DNA-free RNA was re-suspended 
in 200/il of RNase-free water. Intact RNA was confirmed by 
observation of sharp rRNA bands in native agarose gel elec- 
trophoresis. 



RT-qPCR 

Reverse transcription (RT) was performed using Ay\ of the 
DNA-free RNA (100-200ng) and Maxima reverse transcrip- 
tion kit (K1641, Fermentas), used according to the manufac- 
turer's instructions. Random hexamer primers were used for 
priming the RT reaction. At the end of the RT reaction, 
the 20/il reaction was diluted 100-fold and 10/il of this di- 
luted sample was used for qPGR in the next step. qPGR was 
performed using Maxima SYBR-Creen qPCR kit (K0221, Fer- 
mentas) and manufacturer's instructions were followed. qPCR 
was performed in triplicates for each RT reaction and appro- 
priate negative RT controls were used to confirm the absence 
of DNA contamination. gapA mRNA was used as internal ref- 
erence to normalize all other mRNA levels. Standard curves 
with 6 serial dilutions were used to optimize reaction condi- 
tions and ensure amplification efficiency of between 90-100% 
for the yfp and gapA amplicons. AACt method was used to 
obtain the change in mRNA levels due to amino acid limita- 
tion. The qPCR primer sequences are given in Table S8. 



Western blotting 

Fresh colonies were used to inoculate overnight cultures. 
These overnight cultures were then diluted 1:100 into 1ml 
of rich-defined medium with all 20 amino acids (see section 
on growth and fluorescence measurements for media composi- 
tion). After approximately 3.5 hours of growth at 30° C when 
O-Deoo was ~0.4, cells were spun down at 9000g for 1 min, 
and then re-suspended in 1ml of rich-defined medium with- 
out the amino acid whose limitation was to be induced. This 
re-suspended culture was then split into two equal aliquots. 
The limiting amino acid was added to one aliquot (as a rich- 
medium control) while the other aliquot did not have the 
limiting amino acid. The re-suspended medium also con- 
tained 200ng/ml of anhydro-tetracycline in order to induce 
the pitetOl promoter that controls the 3XFLAG-7//p vari- 
ants. After growth at 30° C for 60 min, cells were spun down 
at 12000g, 1 min and re-suspended in 40-400^1 of CellLytic B 
buffer (Sigma, B7435). The buffer volume used was propor- 
tional to the ODqoo measured at the time of harvesting the 
culture. The lysate was stored at -80°C. lO/xl of the lysate 
was mixed with 2X Laemmli Buffer (Biorad) and then loaded 
onto each lane of a pre-cast polyacrylamide gel (Biorad) and 
SDS-PAGE was carried out at lOOV for 120 min. Proteins 
were transferred to a nitrocellulose membrane by semi-dry 
blotting at 180mA for 60 min. The membrane was blocked 
in 2% skim-milk-TBST overnight, and then incubated with 
a 1:2000 dilution of an anti-FLAG antibody (F3165, Sigma) 
in 10ml of 2% skim-milk-TBST with shaking at room tem- 
perature for 90 min. After washing 4 times for 5 min with 
TBST, the membrane was incubated with 1:2000 dilution of a 
secondary HRP-conjugated antibody (7076, Cell Signaling) in 
15ml of 2% skim-milk-TBST with shaking at room tempera- 
ture for 60 min. After washing 4 times for 5 min with TBST, 
the membrane was treated with an HRP substrate (L00354, 
Genscript) for 5 min and exposed for 30s to a luminescence 
imager. 



12 



Subramaniam ct dl. 



Analyses 

Matlab R2009b (Mathworks) was used for all analyses unless 
otherwise mentioned. All correlations and P- values reported 
in this work were calculated using the Matlab command 'corr' 
with the 'Type' option set to either 'Spearman' or 'Pearson' 
as appropriate. 

Growth and fluorescence analysis. Background absorbance 
and fluorescence values (obtained from wells containing only 
growth media) were subtracted from the measured time series 
for each well. An exponential curve was fitted to the amino 
acid-rich growth regime for all data points located at least 
50 min before the onset of amino acid limitation. A straight 
line was fitted to the amino acid-limited growth regime for all 
data points located at least 50 min after the onset time. These 
fits were performed using the Matlab command 'fit', and the 
in-built library options 'Expl ' and 'Polyl ' respectively. To 
automatically identify the onset time, the intersection point 
between the two fitted curves was designated as the onset time 
of amino acid limitation. This inferred onset time coincided 
with the onset time identified through visual inspection of the 
growth curves. 

To minimize noise in calculated protein synthesis rates, an 
exponential curve was fitted to the amino acid-rich regime of 
the fluorescence time-series and a straight line was fitted to 
the amino acid-limited regime of the fluorescence time-series. 
These fits were performed using the Matlab command fit, and 
the in-built library options 'Expl ' and 'Polyl ' respectively. 
Protein synthesis rate, S was calculated as 



Protein synthesis rate S 



d{Fl 



uorescence 



Absorbance 



d(timc) 



[81 



First, the above formula was evaluated at the onset time of 
amino acid limitation using the exponential fits for absorbance 
and fiuorescence data in the amino acid rich growth regime. 
Next, the same formula was evaluated at the onset time us- 
ing the linear fits for absorbance and fiuorescence data in the 
amino acid limited growth regime. These two values corre- 
spond to the protein synthesis rates reported for the amino 
acid rich and amino acid limited growth regimes (such as the 
data in Fig. ID). The protein synthesis rates were normalized 
within each synonymous codon family and for each growth 
condition. Robustness of protein synthesis to amino acid lim- 
itation was calculated as the ratio of normalized protein syn- 
thesis rates between the amino acid rich and amino acid lim- 
ited growth regimes. 

In the case of the experiment without methyl ester (Fig. 
S15 and Fig. S16), the onset time of amino acid limited growth 
was determined exactly as above. Then starvation robustness 
was calculated as the normalized ratio of total fiuorescence 
increase after the onset of amino acid limited growth to the 
fluorescence increase before this onset. Total fluorescence in- 
crease rather than protein synthesis rate was used for this 
analysis since protein synthesis rates decreased continuously 
to zero after the onset of amino acid limited growth in the 
absence of methyl ester analogs. 

Calculation of CRI. CRI for a protein coding sequence corre- 
sponding to a limiting amino acid was calculated by multiply- 
ing the Wi values for codons cognate to the limiting amino 
acid in that sequence. Wi values in Fig. 3B were calcu- 
lated using the robustness of protein synthesis of the cor- 
responding yfp variants during cognate amino acid limita- 
tion (Fig. ID). Based on the non-cognate amino acid lim- 
itation experiment (Fig. S2), the Wi values for all codons 



other than those cognate to the limiting amino acid are 
set to be equal to 1. For illustration, we demonstrate 
the calculation of Wi for the six Leu codons below. The 
exact same procedure was followed for other synonymous 
codon families. Taking logiWi = Wi for each codon, and 
Zo(;2 (robustness during amino acid limited growth) = SR for 
each yfp variant, 

7 X WCTA -H5 X WCTG = SRyfp,CTA [9] 

7 XWCTC + 15 XWCTG = SRyfp,CTC [10] 

22 X WcTG = SRyfp,CTG [11] 

7 XWCTT + 15 xWcTG = SRyfp^CTT [12] 

7 xWtTA + 15 xWcTG = SRyfp^TTA [13] 

7 XWtTG + 15 xWcTG = SRyfp^TTG [14] 

The multiplicative factors on the LHS in front of Wi corre- 
spond to the frequency of the Leu codon i in the corresponding 
Leu variant of yfp (see Fig. lA). The RHS is the measured 
(log2) robustness of protein synthesis from the corresponding 
yfp variant during Leu limitation (see Fig. ID). These equa- 
tions were solved simultaneously to determine the Wi value 
for each Leu codon. Revised Wi values (Table S7) based on 
yfp measurements in the presence of ''^'^Leu2 tRNA (Fig. 2) 
were used for calculation of Leu CRI in the case of '^'^^Leu2 
tRNA co-expression with E. coli ORFs (Fig. 3D). 

Leu and Arg CRI for E. Coli ORFs. 4300 E. coU ORF se- 
quences were parsed out from the MG1655 genome sequence 
(NCBI website. Accession number: NC_000913, downloaded 
on 14th Apr 2011). For each of these 4300 E. coU ORFs, 
Leu or Arg CRI was calculated by multiplying the Wi values 
for either all Leu or all Arg codons respectively in the ORF 
sequence. For the 63 synonymous variants of 13 ORFs (Fig. 
SIO), Leu CRI values were calculated using the same proce- 
dure as above after accounting for the synonymous mutations. 
For the 21 ORFs co-expressed with Leu2 tRNA (Fig. 3D), re- 
vised Wi values were flrst calculated using the method outlined 
in the previous section (Table S7), and using measurements on 
the 6 Leu variants of yfp complemented with ''*'^Leu2 tRNA 
(3rd column in Fig. 2A). These revised Wi values were then 
used to calculate Leu CRI under tRNA co-expression for the 
21 tRNA co-expressed ORFs applying the same procedure as 
for the non co-expressed case. 

Z-score for CRI. To quantify the deviation in CRI from its ex- 
pected value for each of the 4300 ORFs in the E. coli genome, 
1000 random coding sequences were generated for each ORF. 
Each random version preserved the original amino acid se- 
quence, but the codons for a single amino acid were sampled 
randomly from a multinomial distribution based on the aver- 
age frequency of codons for that amino acid in the genome. 
CRI values were calculated for each random version of the 
gene, and a distribution of CRI values was generated from 
the 1000 random trials. The average, ficRi and standard de- 
viation, acRi of this CRI distribution was used to calculate 
the Z-score for CRI as follows: 



ZcRi ~ 



CRIo. 



bserved — fJ-C RI 
CCRI 



[151 



In the case of the Z-score for leucine shown in Fig. 4B, the 
leucine codon CTG was not randomized in the above calcu- 
lation and only the remaining 5 leucine codons: CTA, CTC, 
CTT, TTA, and TTG were randomized. This step is impor- 
tant since CTG, which is read by an abundant tRNA isoaccep- 
tor, is enriched in highly-expressed genes, and such genes will 
show up falsely as perturbation-robust genes because CTG is 
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also the codon that is most robust to leucine limitation in our 
experiments (see Fig. ID). 

Codon-specific bioinformatic measures. Codon usage in Fig. 
S4 was calculated as the average frequency of each codon 
across the genome of E. coli MG1655 (4300 ORFs). tRNA 
concentrations in Fig. S4 were taken from previous work (see 
Table 2 in 1]). Concentrations of all cognate tRNAs for each 
codon were summed together. The codon-tRNA adaptation 
index in Fig. S4 is taken from literature (see Table S2 in ['>'>]). 
The tAI value for the CCA codon was revised from the unreal- 
istically low value of 0.00005 to 0.1333 as explained previously 
[5(i] . For inferring codon elongation rates from charged tRNA 
fractions (Fig. S5), we used the formula for codon elongation 
rate from [!)]: 



Vk Zjitiailtik 

where Vk is the elongation rate of codon k, tq is the codon- 
independent elongation time across any codon, ti is the con- 
centration of tRNA isoacceptor i that is cognate to codon 
k, kik is the second-order rate constant for binding of the 
ternary complex containing the charged isoacceptor i to the 
ribosome at codon k, and Qi is the charged fraction of isoac- 
ceptor i. We calculated the codon elongation rates during 
amino acid limitation using the measured charged fractions 
from [10]. For amino acid rich conditions, we set the charged 
fraction to be equal to unity. We used ro = 0.05s~^ , and 
kik = 2 X 10^M~^s~^ similar to [!)]. The ratio of codon elon- 
gation rates was then normalized within each codon family by 
the maximum value within that family. 

ORF-specific bioinformatic measures. Codon Adaptation In- 
dex was calculated for each E. coli ORF using the method 
in [18]. This calculation was implemented using the Codon- 
Adaptationlndex class in the CodonUsage module of BioPy- 
thon (version 1.58). tRNA Adaptation Index was calcu- 
lated for each E. coli ORF using the method in [57]. This 
calculation was implemented using the codonR package ( 
http : //people . cryst . bbk. ac .uk/-f dosrOl/tAI/ index .html, 
downloaded on 3rd Sep 2011). mRNA folding energy was cal- 
culated for the first 37nt of each E. coli ORF together with 
the 5 upstream nucleotides (GTACC) in the pZSll plasmid 
backbone. Calculation was implemented using the hybrid-ss- 
mm command in UNAFold software v3.8 [58] with default 
parameter values for reaction conditions (NA = RNA, T = 
37, [Na-f] = 1, [Mg++] = 0, maxloop = 30). 

Supplementary notes 



Use of methyl esters in amino acid limitation experiments. 

Before we settled on the methyl ester analog-based experi- 
ments, we tested two other amino acid limitation assays that 
are commonly used in the literature. The first assay is a spin 
— > wash — > resuspend in amino acid-|- / amino acid- medium 
[10]. We did not pursue this assay for most experiments since 
it is logistically difficult to perform this assay when working 
simultaneously with more than a dozen strains. However, we 
used this assay for the Western blotting and RT-qPCR mea- 
surements on a few strains (Figs. SI, S14 and 4B). 

The second assay involves starting with a low initial con- 
centration of an amino acid and letting the bacterial cultures 
exhaust the amino acid in the medium through exponential 
growth [5!)]. The bacteria then enter the amino acid limited 
regime in mid-log phase without any intervention from the ex- 
perimenter. However in the absence of exogenous sources of 
amino acid in the amino acid limited regime, protein synthe- 
sis occurs only transiently for less than an hour under these 
conditions and YFP synthesis rates from all yfp variants drop 
below measurable levels at the end of this time period (Fig. 
S15). More importantly, there is no extended steady state 
during which differential protein synthesis rates can be mea- 
sured accurately. Nevertheless, we have confirmed that the 
measurements with and without methyl esters give qualita- 
tively similar results (Fig. S16). In addition. Western blotting 
done in the absence of methyl ester reproduced the heirarchy 
in protein levels between synonymous variants of YFP during 
amino acid limitation (Fig. SI). 

In contrast to the assay without methyl ester, presence of 
methyl ester analogs in the growth medium results in a quasi- 
steady state of amino acid limited growth due to hydrolysis 
of the ester, during which differential YFP expression can be 
measured easily (Fig. S15). Such partial amino acid limited 
growth is also likely to be the relevant scenario when pro- 
totrophic strains run out of amino acids in their growth me- 
dia and have a limited supply of amino acids through protein 
degradation or partially up-regulated biosynthesis pathways. 

Effect of varying the initial concentrations of amino acids 
and methyl esters. Increasing the initial concentration of the 
amino acid or its methyl ester results in a higher cell density 
for the onset of amino acid limitation, and when the corre- 
sponding concentrations are decreased, this onset happens at 
a lower cell density. Importantly, the observed differential ro- 
bustness of protein synthesis (such as the data shown in Fig. 
ID) is qualitatively the same upon 2-fold changes to the initial 
concentration of either the amino acid or its methyl ester. As 
an extreme example, see Figs. S15 and S16 for comparison 
between the cases with and without methyl ester analog in 
the growth medium. 
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Fig. SI: Expression level of yfp variants quantified through Western blotting 

Modified versions of 29 yfp variants (Fig. lA) were created that had a 3X-FLAG tag at the 5' end. These yfp variants were 
transformed into the respective E. coli auxotrophs in which YFP synthesis was repressed by the TetR protein. Cells were 
harvested at an ODsoo of 0.4 and re-suspended in medium with or without the corresponding amino acid. Expression of YFP 
was induced using 200ng/ml anhydrotetracyline, and cells were harvested after 60 min. For each set of yfp variants under a 
specific growth condition, the same amount of total protein (as measured by ODaoo before cell lysis) was used for Western 
blotting. 
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Fig. S2: Effect of a single synonymous mutation on YFP syntfiesis rate 

22 variants of yfp were synthesized, each of wliicli had a single CTA codon at one of the 22 leucine codon locations along yfp. 
The remaining leucine codons in each variant were the perturbation-robust CTG codon. The 'control' yfp variant did not have 
any CTA codon. Vertical axis refers to the YFP synthesis rate from the 22 variants normalized by that of the control variant, 
either during leucine limitation (top panel) or during leucine-rich growth (bottom panel). Horizontal axis indicates the location 
of the CTA codon along each yfp variant (ATG start codon = 1). Error bars show standard error over three replicate cultures. 
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Fig. S3: YFP synthesis rate during limitation for a non-cognate amino acid 

Leucine and arginine variants of yfp were expressed in an E. coli strain, CP78, that is auxotrophic for both leucine and arginine. 
Response of the 6 Leu variants to Arg limitation is determined by the Arg codons in yfp (CGT and CGC) that are common 
across all 6 Leu variants. Reciprocally, the response of the 6 Arg variants to Leu limitation is determined by the Leu codon 
that is common to the Arg variants of yfp (CTG). YFP synthesis rates are defined as in Fig. ID. Error bars show standard 
error over three replicate cultures. 
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Fig. S4: Comparison of synthesis rate robustness witli codon usage and tRNA concentration 

(A) Codon usage was calculated as the average frequency of each codon across all protein coding sequences in E. coli. (B) 
tRNA concentration for each codon was calculated as the sum of tRNA concentrations for all cognate tRNAs [54]. (C) Since 
tRNAs can differ substantially in their affinity for their cognate codons, we also compared the measured robustness against 
the tRNA adaptation index for each codon [57]. This index accounts for different affinities of synonymous codons for the same 
tRNA isoacceptor. All three measures along the horizontal axes were normalized by the maximum value within each codon 
family. Robustness to amino acid limitation was quantified as the ratio of normalized YFP synthesis rates between amino acid 
limited and amino acid rich growth phases. Error bars represent standard error over three replicate cultures. The data points 
that are not visible for a few codons overlap at the top right-hand corner of each plot. 
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Fig. S5: Comparison of synthesis rate robustness with charged tRNA fraction 

To compare YFP synthesis rates with charged tRNA fractions, the elongation rates for leucine and arginine codons were inferred 
from the measured charged fraction of leucine and arginine tRNA isoacceptors [10] (see section on codon specific bioinformatic 
measures). Previously assigned codon-tRNA assignments and kinetic parameters were used [9]. Note that charged tRNA 
fractions cannot be directly compared with synthesis rates of yfp variants due to overlapping and multiple codon assignments 
for several tRNA isoacceptors. Robustness to amino acid limitation was quantified as the ratio of normalized YFP synthesis 
rates between amino acid limited and amino acid rich growth phases. Error bars represent standard error over three replicate 
cultures. Relative codon elongation rate is the ratio of codon elongation rates between amino acid starved and amino acid rich 
growth regimes, normalized to the maximum value within each synonymous codon family. 
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Fig. S6: Miscoding of a single arginine residue in YFP causes loss of fluorescence 

To test whether mistranslation of arginine residues can underlie the high residual fluorescence of Arg yfp variant-Arg tRNA 
pairs (AGA: argS, AGG: arg3, and CGG:arg4, argS in Fig. 2B), three YFP mutants were created that had one of three single 
point mutations at Arg96: R96H, R96K, and R96Q. The mutant and the 'wild-type' YFP proteins were expressed from a 
pUC18 high-copy vector. Each of the three mutations at Arg96 to a chemically similar amino acid (H, K or Q) decreased YFP 
fluorescence to background level (that of an empty pUC18 vector). Error bars denote standard deviation over five biological 
replicates. 
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Fig. S7: Histogram of CRI 

Green and purple data markers correspond to the Leu and Arg CRI values for 4300 ORFs in E. call's genome. Blue and red data 
markers correspond respectively to Leu and Arg CRI values for the E. coli OKF-yfp fusions that were used to experimentally 
validate CRI. 



18 



Subramaniam ct Cll. 




Arg codon robustness index 

Fig. S8: Correlation of Arg CRI with measured robustness of 56 E. coli ORF-yfp fusions 

The yfp sequence used for this experiment had the AGA codon at all Arg codon locations of yfp since AGA has the highest 
Wi value among arginine codons (see Fig. 3B). Correlation is reported as squared Spearman rank correlation. Error bars show 
standard error over three replicate cultures. Robustness to amino acid limitation was quantified as the ratio of normalized YFP 
synthesis rates between amino acid limited and amino acid rich growth phases. 
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Fig. S9: Correlation of protein synthesis rate during amino acid rich growth with measures of translation efficiency 

Protein synthesis rates from 92 E. coli OKF-yfp fusions during Leu-rich growth showed only a weak correlation with measures 
of codon adaptation, tRNA adaptation and 5' folding energy of mRNA. Folding energy was calculated from -5 to -f 37 nt of the 
ATG codon. Godon adaptation index (CAI) and tRNA adaptation index (tAI) were calculated using Biopython and codonR 
packages. Correlations are reported as squared Spearman rank-correlation coefficient. 



Subramaniam et dl. 



19 



variantl 
variant2 

CTC^ 



wild-type ORF 
— TTA / TTG 



DC 
O 



8 



Uj 




2-5 1 - 

Leu codon 
robustness index 

(Prediction) 



1 
Robustness to 
Leu-limitation 

(Measurement) 



Fig. SIO: CRI predicts the change in robustness during amino acid limitation due to synonymous mutations 

Sixty three synonymous variants of 13 OKF-yfp fusions were constructed by mutating wild-type TTG or TTA codons in 
the ORF sequence to the codon CTC that causes sensitive protein synthesis rate under leucine limitation. The number of 
mutations was between 1 and 6 and the location of these mutations was random. 59 of the 63 variants displayed a decrease 
in their robustness during Leu hmitation (dashed lines) that was predicted by Leu CRI (sohd hues). In addition, magnitude 
of the changes in robustness during Leu limitation were positively correlated with magnitude of the changes in Leu CRI {r^ 
= 0.19, P = 10"*). Filled circles indicate values for ORF iid-typc and open circles indicate values for ORFvariant. Different 
open circles within a single polygon correspond to distinct ORF variants for the same wild-type ORF. Robustness to amino 
acid limitation was quantified as the ratio of normalized YFP synthesis rates between amino acid limited and amino acid rich 
growth phases. Error bars show standard error over three replicate cultures. Most error bars are smaller than data markers. 
DNA sequences for variants are provided in gene_sequences.fasta supplementary file. 
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Fig. Sll: Effect of synonymous mutations and tRNA co-expression on synthesis rate from E. coli ORF-yfp fusions 

We analyzed the change in protein synthesis rates from the 21 OKF-yfp fusions co-expressed with *^^*'Leu2 (Fig. 3D) and 
the 63 different OKF-yfp variants with synonymous mutations (Fig. SIO). Several of the *^*'^Leu2-coexpressed as well as the 
synonymously-mutated ORF-yfp variants (84 total variants) had significantly altered protein synthesis rates compared to their 
non-tRNA co-expressed or non-mutated counterparts (referred as wild type) during leucine limited growth (green histogram, 
median fold-change in protein synthesis rates = 2.37). By comparison, most of the 84 variants had similar protein synthesis 
rates to their wild-type counterparts during leucine rich growth (grey histogram, median fold-change in protein synthesis rates 
— 1.12). Protein synthesis rates were defined as in Fig. ID. 
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Fig. S12: Effect of synonymous mutations on synthesis rate from amino acid biosynthesis genes 

(A) Synthesis rates from leuA, leuC, leuD and carA, argA, argF -yfp fusions encoding either wild-type or mutant ORF sequences 
during amino acid rich and amino acid limited growth. The synthesis rates were normalized for each pair of wild-type and 
mutant OKF-yfp fusions, and also for each growth condition. Error bars show standard error over six replicate cultures. (B) 
Position and identity of synonymous mutations in wild- type and mutant sequences used for the experiment in (A). The black 
vertical bars correspond to the non-mutated leucine codons in the case of leuA, leuC and leuD, and to the non-mutated arginine 
codons in the case of car A, argA and argF. 
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Fig. S13: Correlation of protein synthesis rate with ribosome occupancy at serine codons during serine-limited growth 

Protein synthesis rate of serine synonymous variants of yfp during serine limitation (same data as in Fig. ID, third panel) 
is negatively correlated with genome-wide ribosome occupancy at serine codons during serine-limited growth of E. coli. The 
increased occupancy at perturbation-sensitive serine codons is consistent with selective ribosome pausing at these codons. 
Ribosome occupancy data was taken from a recent ribosome profiling experiment in E. coli [fiO]. 
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Fig. S14: Change in mRNA level of yfp variants in response to cognate amino acid limitation 

We measured the change in mRNA levels of different yjp variants in response to amino acid limitation. Total RNA was 
extracted either during exponential amino acid rich growth or 60 min after amino acid limited growth in the presence of the 
amino acid methyl ester. mRNA levels were quantified by RT-qPCR relative to gapA mRNA. Error bars show standard error 
of triplicate qPCR measurements. Synthesis rate robustness to amino acid limitation was quantified as the ratio of normalized 
YFP synthesis rates between amino acid limited and amino acid rich growth phases. 
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Fig. S15: Raw absorbance and fluorescence curves with or without methyl ester analog of amino acids 

Growth and fluorescence curves for two yfp variants corresponding to CTA and CTG codons are shown here as representative 
examples for amino acid hmited growth in the presence or absence of methyl ester analogs in the growth medium. Absorbance 
as measured using spectrometry is proportional to cell density. Presence of methyl ester analogs caused an increase in the time 
and cell density at which amino acid limited growth began. More importantly, inefficient metabolism of methyl ester analogs 
resulted in a slow but steady growth in amino acid limited regime. This residual growth ensured that YFP synthesis continued 
robustly from the CTG yfp variant under these conditions. By contrast, in the absence of methyl esters in the growth medium, 
YFP synthesis from all yfp variants eventually dropped to zero. 
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Fig. S16: Synthesis rate robustness with or without methyl ester analog of amino acids 

Robustness to amino acid limitation in the absence of methyl ester analogs was calculated as the ratio of fluorescence change 
between the amino acid limited growth phase and amino acid rich growth phase. This ratio was further normalized by the 
maximum value within each codon family. Robustness to amino acid limitation in the presence of methyl ester analogs was 
quantified as the ratio of normalized YFP synthesis rates between amino acid limited and amino acid rich growth phases. Error 
bars show standard error over three replicate cultures. 
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Fig. S17: Plasmid map of expression vector for yfp and E. coli ORF-yfp fusions 

A specific plasmid construct with yfpO is shown here. In the case of ORF-y/p fusions, yfp was fused in-frame to the 3'-end 
of the ORF with a GGSGGS hexa-peptide linker sequence that encoded a BamHI restriction site and the resulting coding 
sequence of the fusion protein was cloned between the Kpnl and Hindlll restriction sites in the above vector. 
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Fig. S18: Plasmid map of expression vector for tRNA genes 

A specific construct encoding an Arg tRNA is shown here. 
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Fig. S19: Reproducibility of measurements between biological replicates of 92 E. coli ORF-yfp fusions 

Two different colonies were picked after cloning the 92 E. coli OKF-yfp fusions and the same leucine limitation assay that was 
used for the data in Fig. 3C was performed on these two biological replicates on two different days. None of the clones for 
replicate 2 were sequence-verified and hence the few outliers seen above could be the result of errors in the cloned sequences. 
The data reported in Fig. 3C is from replicate 1 for which about 40 constructs were sequence-verified. Robustness to Leu 
limitation was calculated as the ratio of normalized YFP synthesis rates between Leu limited and Leu rich growth phases. 



Table SI: 92 E. coli ORF-yfp fusions used for Leu CRI validation 

Genes are arranged by increasing values of Leu CRI. SLeu-rich and SLeu-Umited refer to respective protein synthesis rates (a.u. 
per sec per cell). Robustness refers to the ratio between the two protein synthesis rates after normalization by the corresponding 
value for the CTG variant of yfp (which is the yfp tag in these ORF-y/p fusions). ± refers to standard error of measurement. 

Number Gene SLeu-rich SLeu-Umited Robustness log2(Leu CRI) Gene product 



1 


polB 


10.4 


± 0.3 


0.6 


± 


0.6 


0.063 


± 


0.062 


-19.45 


DNA polymerase II 


2 


thiH 


69.7 


± 1.1 


-0.7 


± 


0.8 


-0.012 


± 


0.013 


-14.52 


tyrosine lyase, involved in 
























thiamin-thiazolemoiety synthe- 
sis 


3 


aat 


26.7 


± 0.8 


1.3 


± 


0.3 


0.054 


± 


0.015 


-12.52 


leucyl / phenylalanyl-tRN A- 
























protein transferase 


4 


gdhA 


20.4 


± 0.3 


-0.2 


± 


0.2 


-0.013 


± 


0.009 


-11.62 


glutamate dehydrogenase, 
























NADP-specific 


5 


serC 


91.3 


± 0.5 


-1.7 


± 


0.6 


-0.020 


± 


0.008 


-11.6 


3- 
























phosphoserine / phosphohydroxy 
























threonine aminotransferase 


6 


gpsA 


111.8 


± 3.3 


1.9 


± 


1.9 


0.020 


± 


0.018 


-11.38 


glycerol-3-phosphate dehydro- 
























genase (NAD-I-) 


7 


mlrA 


44.1 


± 1.5 


-0.7 


± 


0.8 


-0.017 


± 


0.019 


-11.09 


DNA-binding transcriptional 
























regulator 


8 


ybeU 


41.5 


± 7.6 


0.5 


± 


0.3 


0.017 


± 


0.013 


-11 


conserved protein, DUF1266 
























family 


9 


argS 


67.6 


± 1.9 


1.0 


± 


0.6 


0.017 


± 


0.010 


-10.93 


arginyl-tRNA synthetase 


10 


ilvD 


59.2 


± 5.6 


3.4 


± 


0.9 


0.064 


± 


0.017 


-10.39 


dihydroxyacid dehydratase 


11 


UgpC 


88.5 


± 2.7 


1.9 


± 


0.9 


0.023 


± 


0.011 


-10.33 


glycerol-3-phosphate trans- 
























porter subunit 


12 


argI 


44.5 


± 4.5 


5.3 


± 


2.6 


0.138 


± 


0.067 


-10.24 


ornithine carbamoyltransferase 
1 


13 


ttdA 


86.1 


± 0.2 


-0.0 


± 


0.5 


-0.000 


± 


0.006 


-10.06 


L-tartrate dehydratase, alpha 
























subunit 


14 


leuS 


143.4 ± 3.7 


-2.5 


± 


0.5 


-0.019 


± 


0.003 


-9.46 


leucyl-tRNA synthetase 


15 


hisB 


119.8 


± 3.7 


-0.3 


± 


3.0 


-0.002 


± 


0.027 


-9.31 


fusedhistidinol- 



phosphatase / imidazoleglycerol- 
phosphatedehydratase 
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Table SI (contd.): 92 E. coli ORF-yfp fusions used for Leu CRI validation 



Number Gene SLeu-rich SLeu-Hmited Robustness log2(Leu CRI) Gene product 



16 


rpoA 


71.8 ± 7.5 


0.0 


17 




173 2 ± 8 8 


-5.2 


18 


rob 


24.9 ± 0.7 


0.5 


19 


agaS 


73.0 ± 3.3 


1.3 


20 


lysS 


111.8 ± 4.4 


3.7 


21 


sdaB 


121.0 ± 5.7 


36.4 


22 


purA 


35.5 ± 2.2 


-1.2 


23 


rpoD 


77.5 ± 2.2 


2.3 


24 


melR 


44.5 ± 3.3 


15.1 


25 


kefF 


99.2 ± 3.0 


3.3 


26 


uhpA 


22.4 ± 1.5 


0.4 


27 


ydcN 


41.6 ± 2.7 


1.9 


28 


aspS 


62.3 ± 4.2 


10.9 


29 


relB 


87.9 ± 2.3 


0.8 


30 


guaA 


122.6 ± 2.3 


6.9 


31 


phnM 


203.7 ± 9.5 


20.9 


32 


tdcD 


4.5 ± 0.4 


7.1 


33 


ompR 


114.2 ± 18.7 


-1.4 


34 


phnL 


75.7 ± 5.3 


5.7 


35 


purH 


29.6 ± 1.4 


10.2 


36 


argG 


2.5 ± 0.2 


3.0 


37 


rimM 


107.5 ± 3.7 


13.3 


38 


ubiC 


93 6 + 3 9 


7.5 


39 


leuL 


48.9 ± 2.3 


0.7 


40 


asnS 


9.0 ± 0.4 


7.2 


41 


ribB 


107.6 ± 7.8 


12.2 


42 


smpB 


15.8 ± 0.8 


1.5 


43 


guaB 


31.0 ± 0.3 


5.2 


44 


proC 


69.2 ± 3.1 


9.4 


45 


pth 


89.0 ± 2.6 


25.2 


46 


ivbL 


1.8 ± 0.3 


1.6 


47 


chbR 


25.6 ± 1.0 


24.4 


48 


leuA 


78.1 ± 2.2 


53.6 


49 


argD 


85.7 ± 4.4 


51.0 


50 


yfcN 


98.8 ± 2.1 


24.2 



± 0.6 


-0.001 


± 


0.009 


-9, 


,05 




-0.027 




n fifii 


fi 

o. 


■ OKJ 


± 0.2 


0.024 


± 


0.007 


-7, 


,94 


± 0.4 


0.019 


± 


0.005 


-7, 


,48 


± 2.1 


0.035 


± 


0.020 


-7, 


,17 


± 4.1 


0.333 




0.044 


7 


,14 


± 0.4 


-0.039 


± 


0.015 


-7, 


,12 


± 1.1 


0, 


.032 


± 


0.014 


-6, 


,99 


± 0.4 


0, 


.378 




U.UO-L 


R 
U, 


7fi 


± 0.5 


0, 


.038 


± 


0.007 


-6, 


,52 


± 0.7 


0. 


.019 


± 


0.034 


-6, 


,5 


+ 07 

31 U. ( 


0, 


.050 


-1- 

ZIZ 




« 

-u. 


• O 


± 1.5 


0, 


,191 


± 


0.014 


-6, 


,42 


± 0.9 


0, 


.009 


± 


0.010 


-6, 


,09 


± 2.2 


0, 


.061 


± 


0.019 


-5, 


,81 


± 1.2 


0, 


.113 


± 


0.007 


-5, 


,81 


± 0.2 


1, 


.755 


± 


0.163 


-5, 


,76 


± 2.2 


-0.019 


± 


0.027 


-5, 


,74 


± 1.1 


0, 


.082 


± 


0.013 


-5, 


,51 


± 0.4 


0, 


.383 


± 


0.025 


-5, 


,27 


± 0.2 


1, 


.343 


± 


0.150 


-5, 


,07 


± 0.9 


0, 


.137 


± 


0.014 


-4, 


,95 


± 2.2 


0, 


.087 




0.022 


-4, 


S7 


± 0.5 


0, 


.017 


± 


0.011 


-4, 


,75 


± 0.5 


0, 


,884 


± 


0.071 


-4, 


,65 


± 1.0 


0, 


,127 


± 


0.020 


-4, 


,49 


± 0.4 


0, 


,101 


± 


0.023 


-4, 


,39 


± 0.6 


0, 


,184 


± 


0.022 


-4, 


,21 


± 1.6 


0, 


,151 


± 


0.030 


-3, 


,92 


± 0.9 


0, 


,313 


± 


0.018 


-3, 


,1 


± 0.1 


1, 


,085 


± 


0.295 


-2, 


,89 


± 0.8 


1, 


,054 


± 


0.008 


-2, 


,81 


± 0.8 


0, 


,756 


± 


0.019 


-2, 


,77 


± 1.8 


0, 


,657 


± 


0.031 


-2, 


,72 


± 1.5 


0, 


,270 


± 


0.011 


-2, 


,71 



RNA polymerase, alpha subunit 
DNA-binding response regula- 
tor in two-component regula- 
tory system with BarA 
right oriC-binding transcrip- 
tional activator, AraC family 
tagatose-6-phosphate ke- 
tose/aldose isomerase 
lysine tRNA synthetase, consti- 
tutive 

L-serine deaminase II 
adenylosuccinate synthetase 
RNA polymerase, sigma 70 
(sigma D) factor 
DNA-binding transcriptional 
dual regulator 

flavoprotein subunit for the 
KefC potassium efflux system 
DNA-binding response regula- 
tor in two-component regula- 
tory system wtih UhpB 
predicted DNA-binding tran- 
scriptional regulator 
aspartyl-tRNA synthetase 
Qin prophage; bifunctional 
antitoxin of theRelE-RelB 
toxin-antitoxin system/ tran- 
scriptional repressor 
GMP synthetase (glutamine 
aminotransferase) 
carbon-phosphorus lyase com- 
plex subunit 

propionate kinase/acetate ki- 
nase C, anaerobic 
DNA-binding response regula- 
tor in two-component regula- 
tory system with EnvZ 
carbon-phosphorus lyase com- 
plex subunit 

fused IMPcyclohydrolase / 
phosphoribosyl aminoimidazole 
carboxamide formyltransferase 
argininosuccinate synthetase 
16S rRNA processing protein 
chorismate-pyruvate lyase 
leu operon leader peptide 
asparaginyl tRNA synthetase 
3 ,4-dihydroxy- 2-butanone-4- 
phosphate synthase 
trans-translation protein 
IMP dehydrogenase 
pyrroline-5-carboxylate 
reductase,N AD (P)-binding 
peptidyl-tRNA hydrolase 
ilvB operon leader peptide 
rRepressor, chb operon forN,N'- 
diacetylchitobiose utilization 
2-isopropylmalate synthase 
bifunctional acetylor- 
nithine aminotrans- 
ferase /succiny Idiaminopimelate 
aminotransferase 
conserved protein 
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Table SI (contd.): 92 E. coli ORF-yfp fusions used for Leu CRI validation 



imber 


Gene 


S Leu — rich 


S Leu — limited 


Robustness 


log2(Leu CRI) 


Gene product 


51 




7.3 ± 0.3 


9.9 ± 0.4 


1.503 


± 


0.046 


-2.62 


predicted dioxygenase, LigB 


















TP vn 1 1 ^7 


52 


mdt J 


12.5 ± 1.6 


17.5 ± 0.2 


1.595 


± 


0.204 


-2.52 


TTllll^lH T*11 Cf PTTlllY Q^rQ"i"PTn "("VPTIQ- 




















53 


r. 


sn 6 + 9 5 


31.4 ± 2.0 


0.428 


± 


0.018 


-2.39 


TiInI A — nin n i Ti 1"T'fiTTQPT'iT^"f"ion?il 

A^XlJT. LJliHilllti Li. dilOljl 1 IJ LlUilCLl 


















repressor of the aga regulon 


54 


hisA 


22.0 ± 0.9 


27 6 ± 9 


1.392 


± 


0.087 


-1.86 


N-(5'-pliospho-L-ribosyl- 


















TOTTTl i TTl ITir^l-^— a TTl 1 TT C\— 1 — 
iVJi iiiiiiiiiiLJ 1 CLiiiiiiLJ i 


















( '-Till nsnVinribntivl 1-4- 

\U LJiiLJo i^iiVJi i ULJo y i J T: 


















imidazolecarb oxamide iso- 


















iiiCi CLDC 


55 




72 7 + 21 5 


26.4 ± 1.3 


0.475 


± 


0.130 


-1.65 


TiTnbfiblp ^qRIMa pTidnTiiiplpaQP 

^iLJUCLUiC oo ± Li 1 Ji CiiviLJii LHjiCdoC J 


















CRISP-associatedprotein 


56 


rpoH 


89.2 ± 1.8 


50.5 ± 0.9 


0.623 


± 


0.005 


-1.53 


R,NA polymerase, sigma 32 


















(sigma H) factor 


57 


leuD 


128.4 ± 2.7 


70 ± 6 


0.601 


± 


0.018 


-1.48 


3-isopropylmalate dehydratase 


















small subxmit 


58 


din J 


58 9 ± 2 8 


34.4 ± 1.5 


0.648 


± 


0.050 


-1.48 


antitoxin of YafQ-DinJ toxin- 


















antitoxin system 


59 


nuol 


95 2 + 1 8 


30.6 ± 1.4 


1.361 


± 


0.164 


-1.48 


^1 A T) n ■ 1 1 nini 1 T Ti nTi p riYTrlnTpriiir'- 

i N .r^i_y ii. . li LJiLJ liiiiLJiiC LJA-iLiLJi CLi LLLj 


















tase, chain I 


60 


luxS 


109.5 ± 0.7 


70.3 ± 3.9 


0.706 


± 


0.035 


-1.34 


S-ribosylhomocysteine lyase 


61 


leuC 


104 5 ± 6 


54.0 ± 3.4 


0.573 


± 


0.069 


-1.18 


3-isopropylnialate dehydratase 


















large subunit 


62 


pspA 


69 8 ± 2 6 


51.4 ± 0.6 


0.813 


± 


0.039 


-1.18 


regulatory protein for phage- 




















63 


pyrl 


1 38 9 + 1 3 7 


119 5 + 69 


0.906 


± 


0.091 


-1.15 




















i eg, U.id,ljL)i Y o U. U Lliii L 


64 


btuE 


1 49 1 +46 


1 39 9 + 3 6 


0.982 


± 


0.020 


-0.95 


o'li li" a + n i c\t\ p np"m"V"iri ti qp 


65 


nisrB 


153 6 ± 14 7 


131.7 ± 2.8 


0.958 


± 


0.074 


-0.8 


methionine sulfoxide reductase 
B 


66 


coaD 


94.4 ± 0.3 


96 9 + n 8 


1.122 


± 


0.012 


-0.71 


T^o vi+p+l-i pi n p_-ri n ncjTi n ai"p fiH PTi^r- 
udiiLc Liiciiic uiitjo uiidLc di_idi y 


















ly It ransf erase 


67 


sfsB 


72.2 ± 2.4 


18.9 ± 1.5 


0.291 


± 


0.032 


-0.61 


1) IM A - niTirliTiD' ^""r^Ti QPViT^I'iOTi a 1 

J.-'X'i i \ UlilU.illtl Li-L dilOtjl IVJ LltJlldl 


















i^C''i^\'\Ti^\ CiT riT Tn?ili"nc;p"mpi"firinlitiTn 

OAjiii \ OuIjKJL (J1 llldlUJOClllC LdULJliOlll 


P,R 
uo 


nirD 




22.1 ± 1.5 


0.223 


± 


0.027 


n ^19 


ni + rHp T(::lA^^n\ iPtfid TsJ A F^^P^ W_ 
111 Ui 1 Lt; i Kid LIU Ldot;, IN I -T J 11 


















binding, sniallsubunit 


69 


fdnl 


62.4 ± 1.6 


57.8 ± 3.4 


1.019 


± 


0.034 


-0.43 


form at p ^lpl^v^^^^o■pT^ a qp- 


















N, cytochrome B556 


















(gamma) subunit , nitr ate- 


















inducible 


70 


greA 


172.5 ± 12.4 


90.0 ± 1.2 


0.581 


± 


0.045 


-0.38 


transcript cleavage factor 


71 


hupB 


987 5 + 19 6 


116 5 + 40 


0.450 


± 


0.030 


-0.33 


TTT T T) A-niTiHinD' ■("T"?iTi 'JPTI v\— 

LIKJ ^ J_-/X\ JT. LJllUJ-ilit LI dlioiji. IJJ 


















i"ir\Ti PI 1 T'PO'i 1 1 niriT' ript a qi i r>i i n ii" 

\ji\JLL<XL 1 Cgj LiidLWl , UVj Ldo LI. Lf Ulii L 


72 


p-InF, 


1 44 9 + 9 n 


1 49 n + 9 6 


1.086 


± 


0.035 


-0.33 


"1" n i r>Qii 1 Tfii"p •PiT'^i niHp qiiItii t'^t'^i n q— 


















Tpv^i Qp 1 T n r>rl ati pqp i 

ICi. doC I i. llULJ-dilCOC J 


73 


ogrlv 


95 8 + 9 5 


1 45 9 + 3 3 


1.699 


± 


0.153 


-0.33 


T^oc;ii"i^7"p "rpD'iilfii"0"r r\\ o'TTi'^sri" n 

LHJOlLlVC i. C t Llld LUi. Ul J. tilLJWLll 


















I in apT"^"^/^T^ r^T V^^r\c¥y itptip in^"/^ ^"rlp 
I lliDfcii LlUli Ui JT ZiUgi gclit; illLU Lilt; 


















chromosome ) 


74 


rplD 


20.9 ± 1.5 


41.0 ± 1.6 


2.191 


± 


0.226 


-0.33 


50S ribosomal subunit protein 
L4 


75 


dmsB 


116.2 ± 5.3 


156.3 ± 5.1 


1.486 


± 


0.067 


-0.28 


dimethyl sulfoxide reductase, 


















anaerobic, subunitB 


76 


rpsP 


241.7 ± 3.7 


118.5 ± 0.3 


0.540 


± 


0.007 


-0.19 


30S ribosomal subunit protein 


















S16 


77 


rplX 


150.6 ± 3.2 


83.7 ± 3.5 


0.613 


± 


0.031 


-0.19 


50S ribosomal subunit protein 


















L24 


78 


tpiA 


138.5 ± 11.4 


87.1 ± 2.2 


0.700 


± 


0.047 


-0.19 


triosephosphate isomerase 


79 


gapA 


39.5 ± 2.0 


62.6 ± 1.6 


1.757 


± 


0.135 


-0.19 


glyceraldehyde-3-phosphate de- 


















hydrogenase A 


80 


rpsT 


174.7 ± 12.0 


133.7 ± 4.5 


0.849 


± 


0.050 


-0.14 


30S ribosomal subunit protein 
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Table SI (contd.): 92 E. coli ORF-yfp fusions used for Leu CRI validation 



Number Gene SLeu-rich SLeu-Umited Robustness log2(Leu CRI) Gene product 



81 rpsJ 170.2 ± 37.2 142.3 ± 1.6 1.015 ± 0.219 -0.14 SOS ribosomal subunit protein 

SIO 

82 rplT 120.5 ± 8.5 112.6 ± 3.7 1.039 ± 0.077 -0.14 SOS ribosomal subunit protein 

L20 

83 ahpC 57.0 ± 5.0 64.9 ± 1.7 1.267 ± 0.078 -0.14 alkyl hydroperoxide reductase, 

C22 subunit 

84 rpsK 109.9 ± 9.2 139.0 ± 3.5 1.411 ± 0.113 -0.14 30S ribosomal subunit protein 

Sll 

85 tsf 200.4 ± 2.8 103.8 ± 2.1 0.570 ± 0.008 protein chain elongation factor 

EF-Ts 

86 rpIU 173.3 ± 12.4 92.5 ± 4.9 0.595 ± 0.056 50S ribosomal subunit protein 

L21 

87 rpml 158.7 ± 4.0 104.2 ± 1.4 0.725 ± 0.025 50S ribosomal subunit protein 

L35 

88 yjgF 131.9 ± 5.3 98.1 ± 3.0 0.824 ± 0.057 conserved protein, UPF0131 

family 

89 ppiB 133.6 ± 4.5 114.6 ± 2.7 0.945 ± 0.011 peptidyl-prolyl cis-trans iso- 

merase B (rotamaseB) 

90 yjbJ 134.6 ± 4.4 158.7 ± 1.1 1.300 ± 0.035 conserved protein, UPF0337 

family 

91 rpsi 27.5 ± 2.8 78.6 ± 6.0 3.198 ± 0.324 30S ribosomal subunit protein 

S9 

92 rpsF 31.7 ± 1.6 97.1 ± 0.8 3.392 ± 0.183 30S ribosomal subunit protein 

S6 



Table S2: 56 E. coli ORF-yfp fusions used for Arg CRI validation 

Genes are arranged by increasing values of Arg CRI. SArg-rich and SArg-HmUed refer to respective protein synthesis rates (a.u. 
per sec per cell). Robustness refers to the ratio between the two protein synthesis rates after normalization by the corresponding 
value for the AGA variant of yfp (this AGA variant was also used as the yfp tag in these ORF-j//p fusions). ± refers to standard 
error of measurement. 



Number Gene SArg-ricyi SArg-Umited Robustness log2(Arg CRI) Gene product 



leuS 100.5 ± 2.6 16.1 ± 1.7 0.107 ± 0.012 -10.05 
phoR 16.4 ± 3.2 0.7 ± 0.6 0.027 ± 0.029 -8.52 



rfnG 61.0 ± 3.4 10.8 ± 2.3 0.118 ± 0.022 -8.51 



4 


asnS 


15. 


,5 


± 


1.8 


2.6 ± 0.1 


0, 


,114 


± 


0, 


,014 


-8, 


,38 


5 


guaA 


91, 


,5 


± 


4.2 


11.7 ± 3.7 


0, 


,083 


± 


0, 


,024 


-8, 


,25 


6 


thiH 


53, 


,9 


± 


1.1 


4.4 ± 1.3 


0, 


,054 


± 


0, 


,015 


-7, 


,16 


7 


fruR 


44, 


,7 


± 


1.5 


5.3 ± 0.1 


0, 


,079 


± 


0, 


,004 


-6, 


,54 


8 


argG 


10, 


,0 


± 


0.9 


2.5 ± 1.1 


0, 


,175 


± 


0, 


,088 


-6, 


,49 


9 


gpsA 


81, 


,9 


± 


3.9 


9.1 ± 4.2 


0, 


,073 


± 


0, 


,032 


-5, 


,81 


10 


rpoA 


60, 


,9 


± 


1.6 


10.5 ± 2.1 


0, 


,114 


± 


0, 


,022 


-5, 


,36 


11 


yiiD 


31, 


,2 


± 


1.8 


9.8 ± 2.9 


0, 


,218 


± 


0, 


,078 


-5, 


,05 


12 


agaR 


53, 


,2 


± 


1.8 


2.6 ± 0.8 


0, 


,032 


± 


0, 


,010 


-5, 


,01 


13 


rimK 


17, 


,9 


± 


2.6 


9.8 ± 1.5 


0, 


,393 


± 


0, 


,106 


-4, 


,94 



leucyl-tRNA synthetase 
sensory histidine kinase in two- 
componentregulatory system 
with PhoB 

fused DNA-binding response 
regulator intwo-component reg- 
ulatory system with GlnL: re- 
sponseregulator/sigma54 inter- 
action protein 

asparaginyl tRNA synthetase 
GMP synthetase (glutamine 
aminotransferase) 
tyrosine lyase, involved in 
thiamin-thiazolemoiety synthe- 
sis 

DNA-binding transcriptional 
dual regulator 

argininosuccinate synthetase 
glyceroI-3-phosphate dehydro- 
genase (NAD-I-) 
RNA polymerase, alpha subunit 
predicted acetyltransferase 
DNA-binding transcriptional 
repressor of the agaregulon 
ribosomal protein S6 modifica- 
tion protein 
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Table S2 (contd.): 56 E. coli ORF-yfp fusions used for Arg CRI validation 



Number Gene SArg-rich SArg-Umited Robustness log2(Arg CRI) Gene product 



14 


rpoH 


75.2 


± 2.8 


9.3 ± 2.9 


0.081 


± 


0.023 


-4.76 


RNA polymerase, sigma 32 




















(sigma Hj factor 


15 


melR 


43.7 


± 0.4 


11.1 ± 2.3 


0.170 


± 


0.036 


-4.46 


DNA-bindmg transcriptional 




















dual regulator 


16 


tdcB 


32.1 


± 2.0 


10.2 ± 2.9 


0.208 


± 


0.051 


-3.83 


catabolic threonine dehy- 




















dratase, FLF-dependent 


17 


serC 


68.2 


± 2.3 


28.9 ± 4.2 


0.284 


± 


0.043 


-3.79 


3- 




















phosphoserine /phosphohydroxy 




















threonine aminotransferase 


18 


rnc 


65.7 


± 2.5 


15.4 ± 2.6 


0.157 


± 


0.028 


-3.74 


RNase III 


19 


chbR 


22.3 


± 1.2 


6.2 ± 2.1 


0.190 


± 


0.064 


-3.66 


T> "U'L. £ IVT AT' 

rRepressor, cnb operon forJN,N - 




















diacetylchitobiose utilization 


20 


rsuA 


100.3 


± 5.9 


12.8 ± 2.8 


0.084 


± 


0.015 


-3.65 


16S rRNA U516 pseudouridine 




















synthase 


21 


tauC 


70.1 


± 3.2 


7.4 ± 1.3 


0.071 


± 


0.016 


-3.5 


taurine transporter subunit 


22 


smpB 


55.3 


± 1.3 


17.5 ± 0.8 


0.211 


± 


0.010 


-3.42 


trans-translation protein 


23 


carA 


28.2 


± 1.1 


14.0 ± 0.8 


0.331 


± 


0.005 


-3.4 


carbamoyl phosphate syn- 




















thetase small subunit,glutamine 




















amidotransferase 


24 


ubiC 


70.6 


± 1.8 


4.7 ± 1.9 


0.044 


± 


0.018 


-3.39 


chorismate-pyruvate lyase 


25 


ptn 


71.2 


± 2.4 


20.0 ± 0.8 


0.187 




0.001 


-3.21 


peptidyl-tRiN A hydrofase 


26 


rpsF 


26.6 


± 1.9 


14.7 ± 2.3 


0.366 


± 


0.043 


-3.06 


30a ribosomaf subumt protein 
S6 


27 


gapA 


34.3 


± 0.9 


12.1 ± 1.9 


0.234 


± 


0.030 


-2.86 


gfyceraldehyde-3-phosphate de- 




















hydrogenase A 


28 


yihL 


43.2 


± 3.4 


6.7 ± 1.1 


0.105 


± 


0.022 


-2.85 


predicted DNA-binding tran- 




















scriptional regulator 


29 


allR 


36.0 


± 1.5 


17.1 ± 0.8 


0.318 


± 


0.026 


-2.82 


1~\ AT A T_ • J • J- • J- • 1 

DN A- binding transcriptional 




















repressor for all(allantoin) 




















and gel (glyoxylate) 




















operons;glyoxylate-induced 


30 


yfcN 


74.2 


± 4.8 


38.9 ± 7.3 


0.353 


± 


0.068 


-2.8 


conserved protein 


31 


fdnl 


42.7 


± 1.8 


17.8 ± 3.2 


0.277 


± 


0.042 


-2.65 


formate dehydrogenase- 




















N, cytochrome B556 




















(gamma)subunit, nitrate- 




















inducible 


32 


ruvA 


53.2 


± 4.1 


26.7 ± 1.1 


0.341 


± 


0.040 


-2.59 


component of RuvABC resolva- 




















some, regulatorysubunit 


33 


adiY 


47.5 


± 8.2 


6.1 ± 0.9 


0.086 


± 


0.003 


-2.55 


1~\ AT A I™ • J • J- • J- • 1 

DNA-bmdmg transcriptional 




















activator 


34 


holD 


15.1 


± 3.6 


2.9 ± 1.3 


0.153 


± 


0.066 


-2.43 


T~vAT A 1 TTT " 1 

DInA polymerase 111, psi sub- 




















unit 


35 


dmsB 


68.1 


± 2.2 


54.7 ± 6.5 


0.537 


± 


0.067 


-2.42 


dimethyl sulfoxide reductase, 




















anaerobic, subunitB 


36 


bglJ 


10.5 


± 1.3 


6.6 ± 0.7 


0.420 


± 


0.010 


-2.39 


1~v AT A 1 • 1 ■ 1 ' 1 ' 1 

DN A- binding transcriptionai 




















activator for silentbgl operon. 




















requires the bglJ4 allele to 




















function; LuxRfamily 


37 


yidi 


71.4 


± 1.4 


17.6 ± 0.7 


0.164 




0.005 


-2.32 


OFb-53 (KpLJilj prophage; 




















predicted protein 


38 


argF 


64.9 


± 2.8 


35.7 ± 2.5 


0.371 


± 


0.044 


-2.01 


ornithine carbamoyltransferase 




















2, chain F; CP4-6prophage 


39 


rplU 


102.0 


± 7.4 


20.5 ± 4.2 


0.134 


± 


0.029 


-1.89 


50S ribosomal subunit protein 




















L21 


40 


luxS 


90.3 


± 1.4 


54.9 ± 2.4 


0.406 


± 


0.018 


-1.77 


S-ribosylhomocysteine lyase 


41 


coaD 


70.9 


± 1.5 


45.3 ± 3.4 


0.426 


± 


0.036 


-1.75 


pantetheine-phosphate adeny- 




















lyltransferase 


42 


glnB 


199.0 


± 6.0 


61.9 ± 7.2 


0.207 


± 


0.021 


-1.75 


regulatory protein P-II for glu- 




















taminesynthetase 


43 


uidR 


32.2 


± 1.7 


35.7 ± 1.5 


0.745 


± 


0.064 


-1.5 


DNA-binding transcriptional 



repressor 
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Table S2 (contd.): 56 E. coli ORF-yfp fusions used for Arg CRI validation 



Number Gene SArg-rich SArg-Hmited Robustness log2(Arg CRI) Gene product 



44 


relB 


69.6 ± 4.3 


on o 1 ri c 

29.8 ± 0.5 


0.287 




0.016 


-1.35 


Qin prophage; bifunctional 


















antitoxin ot tneKelli-KelB 


















toxin-antitoxin system/ tran- 


















scriptionalrepressor 


45 


btuE 


106.1 ± 1.4 


52.7 ± 4.0 


0.331 




0.021 


-1.35 


glutathione peroxidase 


46 


argR 


37.3 ± 3.4 


13.7 ± 4.4 


0.235 


± 


0.054 


-1.19 


DN A- binding transcriptional 


















dual regulator, L-arginine- 


















binding 


47 


OgrK 


53.3 ± 1.5 


on o inn 

29.3 ± 0.9 


0.367 




0.013 


-1.14 


positive regulator of P2 growth 


















/ ■ J ■ r r\ 'til 

(insertion ot F2ogr gene into the 


















chromosome) 


48 


hupB 


127.5 ± 3.6 


71.2 ± 1.0 


0.373 


± 


0.008 


-1.13 


HU, DNA-bmdmg transcrip- 


















tional regulator, betasubunit 


49 


ppiB 


95.0 ± 2.7 


57.2 ± 10.5 


0.401 


± 


0.069 


-1.13 


peptidyl-prolyl cis-trans iso- 


















merase B (rotamaseB) 


50 


yjgF 


92.2 ± 4.7 


74.4 ± 3.6 


0.540 


± 


0.032 


-1.13 


conserved protein, UPF0131 


















family 


51 


ke£F 


68.6 ± 2.5 


69.7 ± 7.6 


0.677 


± 


0.058 


-1.11 


flavoprotein subunit for the 


















KefC potassiumefflux system 


52 


yjbJ 


92.0 ± 5.1 


96.6 ± 30.0 


0.730 


± 


0.272 


-1.02 


conserved protein, UPF0337 


















family 


53 


leuL 


36.7 ± 1.4 


31.8 ± 2.2 


0.577 


± 


0.023 


-0.97 


leu operon leader peptide 


54 


rimM 


86.4 ± 1.7 


89.9 ± 3.3 


0.694 


± 


0.029 


-0.82 


16S rRNA processing protein 


55 


ydjo 


9.9 ± 0.9 


20.7 ± 1.6 


1.432 


± 


0.219 


-0.74 


predicted protein 


56 


mdtj 


21.0 ± 0.4 


26.4 ± 1.5 


0.841 


± 


0.060 


-0.41 


multidrug efflux system trans- 



porter 



Table S3: 21 E. coli ORF-yfp fusions co-expressed with '^^°Leu2 tRIMA 

Genes are arranged by increasing values of Leu CRI as calculated for '~'^'-'Leu2 tRNA co-expression. Sl^^u — rich and OLeu — limited 
refer to respective protein synthesis rates (a.u. per sec per cell) under *'^*'Leu2 tRNA co-expression. Robustness refers to 
the ratio between the two protein synthesis rates after normahzation by the corresponding value for the CTG variant of yfp 
(see Fig. ID). ± refers to standard error of measurement. Refer to Table SI for corresponding values without ''^'^Leu2 tRNA 
co-expression. 



Number Gene SLeu-rich SLeu-Umited Robustness log2(Leu CRI) Gene product 



1 


ygiD 


8.2 ± 0.6 


3.8 ± 0.8 


0, 


.515 


± 


0, 


.103 


-5, 


.69 


predicted dioxygenase, LigB 
























family 


2 


chbR 


26.5 ± 2.5 


7.3 ± 1.5 


0, 


.299 


± 


0, 


.043 


-5, 


.68 


rRepressor, chb operon forN,N'- 
























diacetylchitobiose utilization 


3 


yfcN 


92.7 ± 2.6 


12.1 ± 2.9 


0, 


.142 


± 


0, 


.030 


-5, 


.47 


conserved protein 


4 


mdtJ 


12.7 ± 0.8 


5.7 ± 1.7 


0, 


.513 


± 


0, 


.165 


-4, 


.46 


multidrug efflux system trans- 
























porter 


5 


ilvD 


62.5 ± 1.5 


23.6 ± 0.8 


0, 


.417 


± 


0, 


.022 


-4, 


.29 


dihydroxyacid dehydratase 


6 


aspS 


54.1 ± 2.0 


18.5 ± 1.8 


0, 


.379 


± 


0, 


.046 


-4, 


.24 


aspartyl-tRNA synthetase 


7 


lysS 


106.3 ± 4.7 


35.2 ± 0.8 


0, 


.367 


± 


0, 


.025 


-4, 


.22 


lysine tRNA synthetase, consti- 
























tutive 


8 


leuC 


100.7 ± 0.5 


19.0 ± 1.7 


0, 


.208 


± 


0, 


.019 


-4, 


,04 


3-isopropylmalate dehydratase 
























large subunit 


9 


ygbF 


69.4 ± 3.4 


8.1 ± 0.4 


0, 


.129 


± 


0, 


.009 


-3, 


.99 


probable ssRNA endonuclease. 
























CRISP-associated protein 


10 


pspA 


65.1 ± 1.2 


25.9 ± 2.0 


0, 


.437 


± 


0, 


.025 


-3, 


.91 


regulatory protein for phage- 
























shock-protein operon 


11 


guaA 


122.0 ± 2.8 


51.1 ± 1.9 


0, 


.462 


± 


0, 


.018 


-3, 


.27 


GMP synthetase (glutamine 
























aminotransferase) 


12 


btuE 


137.7 ± 4.6 


57.5 ± 3.9 


0, 


.462 


± 


0, 


.044 


-3, 


.23 


glutathione peroxidase 


13 


purA 


73.3 ± 1.8 


34.8 ± 2.7 


0, 


.522 


± 


0, 


.035 


-3, 


,16 


adenylosuccinate synthetase 


14 


ompR 


116.4 ± 1.0 


50.7 ± 3.2 


0, 


.479 


± 


0, 


.027 


-2, 


.86 


DNA-binding response regula- 



tor in two-component regula- 
tory system with EnvZ 
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Table S3 (contd.): 21 E. coli ORF-yfp fusions co-expressed with '^^'^Leul tRNA 



Number Gene SLeu-rich S Leu-Umited Robustness log2(Leu CRI) Gene product 



15 


phnM 


68.2 ± 4.4 


48.4 ± 1.4 


0.786 


± 


0.044 


-2.72 


carbon-phosphorus lyase com- 


















plex subunit 


16 


msrB 


113.4 ± 3.1 


37.9 ± 2.1 


0.370 


± 


0.028 


-2.68 


methionine sulfoxide reductase 


17 


purH 


71.4 ± 1.3 


7.9 ± 2.0 


0.120 


± 


0.028 


-2.61 


B 

fused IMPcyclohydro- 


















lase/phosphoribosyl aminoim- 


















idazole carboxamide formyl- 


















transferase 


18 


coaD 


98.3 ± 7.4 


46.7 ± 1.9 


0.526 


± 


0.018 


-2.58 


pantetheine-phosphate adeny- 


















lyltransferase 


19 


rimM 


114.3 ± 5.9 


51.3 ± 1.1 


0.498 


± 


0.030 


-1.99 


16S rRNA processing protein 


20 


relB 


79.9 ± 3.0 


66.6 ± 3.3 


0.923 


± 


0.077 


-1.73 


Qin prophage; bifunctional 


















antitoxin of theRelE-RelB 


















toxin-antitoxin system/ tran- 


















scriptional repressor 


21 


proC 


78.1 ± 2.0 


36.7 ± 4.4 


0.519 


± 


0.070 


-1.72 


pyrroline-5-carboxylate 



reductase,N AD (P)-binding 



Table S4: 63 synonymous mutants of 13 E. coli ORF-yfp fusions 

SLeu-rich and SLeu-iimited refer to respective protein synthesis rates (a.u. per sec per cell) of mutant ORFs. Robustness refers 
to the ratio between the two protein synthesis rates after normalization by the corresponding value for the CTG variant of yfp 
(see Fig. ID). ± refers to standard error of measurement. Refer to Table SI for corresponding values of wild-type ORFs. The 
DNA sequence of the mutant variants below is provided in the gene_sequences.fasta file. Three of the sequences below did not 
have any mutations compared to the wild-type ORF and were included as internal controls. 



Number Gene SLeu-rich SLeu-Umited Robustness log2(Leu CRI) 



1 


btuE-yfp 


mutant 


1 


137.0 ± 5.3 


50.1 ± 1.1 


0, 


,403 


± 


0, 


,007 


-2.22 


2 


btuE-yfp 


mutant 


4 


131.5 ± 3.5 


20.5 ± 0.5 


0, 


,172 


± 


0, 


,009 


-3.06 


3 


btuE-yfp 


mutant 


5 


138.9 ± 1.8 


34.1 ± 0.1 


0, 


,270 


± 


0, 


,003 


-2.62 


4 


btuE-yfp 


mutant 


7 


123.4 ± 4.9 


48.2 ± 1.8 


0, 


,433 


± 


0, 


,032 


-2.22 


5 


btuE-yfp 


mutant 


8 


130.7 ± 5.2 


28.8 ± 0.9 


0, 


,243 


± 


0, 


,015 


-3.06 


6 


btuE-yfp 


mutant 


9 


128.9 ± 0.9 


40.0 ± 0.6 


0, 


,342 


± 


0, 


,003 


-2.62 


7 


chbR-yfp 


mutant 


1 


29.5 ± 1.7 


17.7 ± 0.3 


0, 


,663 


± 


0, 


,028 


-2.81 


8 


chbR-yfp 


mutant 


2 


26.7 ± 2.2 


15.1 ± 0.4 


0, 


,634 


± 


0, 


,072 


-3.25 


9 


chbR-yfp 


mutant 


3 


31.7 ± 0.9 


11.3 ± 0.6 


0, 


,393 


± 


0, 


,029 


-5.31 


10 


chbR-yfp 


mutant 


4 


36.0 ± 2.5 


9.0 ± 0.3 


0, 


,279 


± 


0, 


,027 


-6.2 


11 


chbR-yfp 


mutant 


5 


28.9 ± 2.4 


10.7 ± 0.2 


0, 


,414 


± 


0, 


,033 


-5.76 


12 


chbR-yfp 


mutant 


6 


32.6 ± 2.7 


5.7 ± 0.5 


0, 


,197 


± 


0, 


,021 


-7.03 


13 


chbR-yfp 


mutant 


7 


32.1 ± 1.2 


7.7 ± 0.5 


0, 


,267 


± 


0, 


,023 


-5.8 


14 


chbR-yfp 


mutant 


9 


31.3 ± 0.9 


7.9 ± 0.5 


0, 


,278 


± 


0, 


,026 


-6.64 


15 


coaD-yfp 


mutant 


1 


92.3 ± 1.3 


39.1 ± 1.3 


0, 


,467 


± 


0, 


,021 


-1.89 


16 


coaD-yfp 


mutant 


2 


91.7 ± 4.1 


55.5 ± 0.8 


0, 


,669 


± 


0, 


,034 


-1.1 


17 


coaD-yfp 


mutant 


3 


94.0 ± 2.4 


47.7 ± 2.3 


0, 


,558 


± 


0, 


,014 


-1.5 


18 


leuA-yfp 


mutant 


6 


70.5 ± 2.1 


11.3 ± 0.3 


0, 


,177 


± 


0, 


,010 


-4.4 


19 


leuA-yfp 


mutant 


8 


77.9 ± 2.2 


6.6 ± 1.4 


0, 


,093 


± 


0, 


,020 


-5.23 


20 


leuA-yfp 


mutant 


9 


70.0 ± 2.1 


8.4 ± 0.2 


0, 


,132 


± 


0, 


,007 


-4.84 


21 


leuB-yfp 


mutant 


7 


273.5 ± 8.7 


6.1 ± 0.9 


0, 


,025 


± 


0, 


,003 


-4.81 


22 


leuB-yfp 


mutant 


8 


273.3 ± 10.6 


9.6 ± 1.0 


0, 


,039 


± 


0, 


,005 


-4.02 


23 


leuB-yfp 


mutant 


9 


274.9 ± 8.1 


8.4 ± 1.7 


0, 


,034 


± 


0, 


,007 


-3.97 


24 


leuC-yfp 


mutant 


5 


86.8 ± 0.6 


8.5 ± 1.7 


0, 


,108 


± 


0, 


,021 


-3.64 


25 


leuC-yfp 


mutant 


7 


95.6 ± 1.9 


5.1 ± 1.7 


0, 


,059 


± 


0, 


,019 


-3.3 


26 


leuC-yfp 


mutant 


8 


100.5 ± 1.7 


4.3 ± 0.9 


0, 


,047 


± 


0, 


,009 


-3.69 


27 


leuC-yfp 


mutant 


9 


80.0 ± 4.9 


14.3 ± 0.8 


0, 


,198 


± 


0, 


,013 


-3.25 


28 


leuD-yfp 


mutant 


2 


123.6 ± 0.9 


18.3 ± 1.2 


0, 


,163 


± 


0, 


,011 


-2.76 


29 


leuD-yfp 


mutant 


3 


123.1 ± 3.7 


39.5 ± 1.0 


0, 


,355 


± 


0, 


,020 


-1.48 


30 


leuD-yfp 


mutant 


4 


122.0 ± 1.0 


32.9 ± 0.3 


0, 


,297 


± 


0, 


,002 


-1.92 


31 


mdtj-yfp 


mutant 


1 


17.3 ± 0.6 


9.1 ± 0.4 


0, 


,580 


± 


0, 


,020 


-4.15 


32 


mdtj-yfp 


mutant 


2 


12.2 ± 0.5 


5.8 ± 0.3 


0. 


,520 


± 


0, 


,008 


-3.31 


33 


mdtJ-yfp 


mutant 


4 


10.0 ± 0.9 


4.3 ± 0.1 


0. 


,475 


± 


0, 


,047 


-4.54 
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Table S4 (contd.): 63 synonymous mutants of 13 E. coli ORF-yfp fusions 



Number 



Gene 



Sl 



S_ 



Leu — limited 



Robustness 



log2(Leu CRI) 



34 


mdtj-yfp mutant 


5 


13.6 ± 2.1 


11.4 ± 0.5 


0.979 


± 


0.190 


-4.15 


35 


mdtj-yfp mutant 


6 


6.4 ± 0.9 


5.1 ± 0.1 


0.903 


± 


0.125 


-3.75 


36 


mdtJ-yfp mutant 


7 


5.9 ± 3.1 


4.3 ± 0.3 


1.251 


± 


0.488 


-4.15 


37 


mdtJ-yfp mutant 


8 


10.2 ± 1.5 


8.3 ± 0.1 


0.932 


± 


0.117 


-4.19 


38 


mdtJ-yfp mutant 


9 


10.3 ± 0.8 


9.6 ± 0.2 


1.049 


± 


0.087 


-3.36 


39 


msrB-yfp mutant 


1 


76.6 ± 1.9 


23.7 ± 0.8 


0.342 


± 


0.014 


-2.48 


40 


msrB-yfp mutant 


2 


65.3 ± 1.6 


20.7 ± 0.8 


0.348 


± 


0.005 


-2.92 


41 


msrB-yfp mutant 


3 


99.1 ± 2.7 


27.8 ± 2.5 


0.310 


± 


0.031 


-2.52 


42 


msrB-yfp mutant 


4 


94.0 ± 1.5 


36.2 ± 0.4 


0.424 


± 


0.005 


-2.08 


43 


pspA-yfp mutant 


1 


68.2 ± 1.9 


16.3 ± 0.6 


0.264 


± 


0.006 


-2.81 


44 


pspA-yfp mutant 


2 


64.2 ± 0.5 


15.0 ± 0.7 


0.257 


± 


0.011 


-2.85 


45 


pspA-yfp mutant 


3 


62.3 ± 1.5 


28.5 ± 0.6 


0.505 


± 


0.016 


-1.97 


46 


pspA-yfp mutant 


5 


73.3 ± 1.3 


14.0 ± 0.9 


0.211 


± 


0.018 


-3.2 


47 


pspA-yfp mutant 


6 


70.1 ± 3.0 


18.3 ± 0.9 


0.288 


± 


0.013 


-2.76 


48 


pspA-yfp mutant 


7 


73.5 ± 1.4 


8.5 ± 0.2 


0.127 


± 


0.002 


-3.64 


49 


pspA-yfp mutant 


8 


73.0 ± 1.2 


12.3 ± 0.8 


0.185 


± 


0.012 


-3.25 


50 


yfcN-yfp mutant 


1 


103.0 ± 3.5 


8.3 ± 1.3 


0.088 


± 


0.012 


-4.82 


51 


yfcN-yfp mutant 


5 


86.0 ± 3.6 


14.9 ± 1.7 


0.192 


± 


0.022 


-3.1 


52 


yfcN-yfp mutant 


6 


90.5 ± 1.6 


5.6 ± 0.7 


0.069 


± 


0.009 


-5.31 


53 


yfcN-yfp mutant 


8 


91.5 ± 4.1 


14.1 ± 1.7 


0.170 


± 


0.017 


-2.71 


54 


yfcN-yfp mutant 


9 


85.5 ± 1.1 


5.2 ± 0.4 


0.067 


± 


0.005 


-4.87 


00 


ygbF-yfp mutant 


o 

Z 


OU.D ± Z.O 


D.Z ± U.o 


U.UoO 


1 

zh 


U.UUD 


-o./o 


56 


ygbF-yfp mutant 


3 


68.5 ± 3.9 


25.4 ± 1.4 


0.413 


± 


0.044 


-1.65 


57 


ygbF-yfp mutant 


5 


79.0 ± 4.1 


25.1 ± 0.2 


0.352 


± 


0.017 


-2.1 


58 


ygbF-yfp mutant 


6 


65.3 ± 2.9 


10.8 ± 0.6 


0.183 


± 


0.019 


-2.44 


59 


ygbF-yfp mutant 


7 


80.5 ± 4.0 


8.3 ± 0.2 


0.113 


± 


0.004 


-3.33 


60 


ygbF-yfp mutant 


8 


70.4 ± 3.1 


10.8 ± 0.9 


0.170 


± 


0.019 


-2.89 


61 


ygiD-yfp mutant 


1 


11.0 ± 1.0 


2.1 ± 0.2 


0.216 


± 


0.009 


-5.17 


62 


ygiD-yfp mutant 


3 


11.3 ± 0.9 


4.2 ± 0.3 


0.417 


± 


0.056 


-4.73 


63 


ygiD-yfp mutant 


4 


10.6 ± 0.5 


5.4 ± 0.6 


0.572 


± 


0.090 


-4.34 


64 


ygiD-yfp mutant 


5 


7.2 ± 1.1 


3.3 ± 0.5 


0.514 


± 


0.041 


-5.17 


65 


ygiD-yfp mutant 


6 


6.8 ± 1.2 


2.7 ± 0.6 


0.476 


± 


0.131 


-4.29 


66 


ygiD-yfp mutant 


7 


12.3 ± 2.5 


2.7 ± 0.4 


0.251 


± 


0.016 


-6.06 



Table S5: List of strains 



Limiting AA Strain designation CGSC number 



Genotype 



Leu, Arg 
Ser 
Pro 
He 
Gin 
Phe 



CP78 
JW2880-1 
JW0233-2 
JW3745-2 
JW3841-1 
JW2580-1 
MG1655 



4695 W3110, argH- leuB- thr- his- thi- 

10234 BW25113, AserA 

8468 BW25113, Aproyl 

10733 BW25113, Ai/uyl 

10775 BW25113, /\glnA 

10048 BW25113, Ap/ieA 

6300 wild-type strain witli Icnown genome sequence 



Table S6: Concentration of amino acids and methyl esters used for amino acid limitation experiments 



Amino Acid 



Amino acid concen- 
tration in overnight 
cultures (/iA/) 



Amino acid concen- 
tration for amino 
acid limitation exper- 
iments (/iA/) 



Amino acid methyl es- 
ter concentration for 
amino acid limitation 
experiments (fiM) 



Catalog number for 
amino acid methyl es- 
ter 



Leu 
Arg 
Ser 
Pro 
He 



800 

800 

10000 

800 

800 



100 
150 
5000 


100 



160 
160 
800 
160 
160 



L1002 (Sigma) 
11030 (Sigma) 
412201 (Sigma) 
287067 (Sigma) 
10522 (VWR) 
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Table S6 (contd.): Concentration of amino acids and methyl esters used for amino acid limitation experiments 



Amino Acid 



Amino acid concen- 
tration in overnight 
cultures (i-lAI) 



Amino acid concen- 
tration for amino 
acid limitation exper- 
iments (fiM) 



Amino acid methyl es- 
ter concentration for 
amino acid limitation 
experiments (fiM) 



Catalog number for 
amino acid methyl es- 
ter 



Gin 
Phc 



800 
800 



400 
50 



400 
50 



68604 (Astatech) 
P17202 (Sigma) 



Table S7: Wi values for Leu codons under '^■^*^Leu2 co-expression 
Codon - log Wi (with *^^'^Leu2 co-expression) 



CTA 
CTC 
CTG 
CTT 
TTA 
TTG 



0.48 
0.05 
0.04 
0.02 
0.26 
0.33 



Table S8: qPCR primer sequences 



Gene 



Forward primer 



Reverse Primer 



yfp TCATGCTGTTTCATGTGATC 

gapA GCTGAAGGCGAAATGAAAGG 

rpoD TGAAGCGAACTTACGTCTGG 

rpoS TCTCAACATACGCAACCTGG 

rpoH TCGTAATTATGCGGGCTATGG 

rpoE CCAGAAGGGAGATCAGAAAGC 

rpoN TGATCCAACTCTCCCAATTCG 

feci ACTACCACAGCTTCCTTAACG 

fliA CGCTATGCTGGATGAACTTCG 

leuA GTCGCTAACTACAACGGTCG 

ilvM GTTTCCACGTCTGCTCAATG 

ilvB TGAGTTTCCGTGTCCAATCC 



AGGGTGATGCTACTTATGGC 
GTACCAGGATACCAGTTTCACG 
AGAACTTGTAACCACGGCG 
AGCTTATGGGACAACTCACG 
CAGTGAACGGCGAAGGAG 
TACCACATCGGGAACATCAC 
TCGTGATTGGCTAACAGATCG 

TTTCGCTGACCATTACCCG 
CTAAACGTTCCGCTACCTCAG 
GCACGCCAGATATTGTTCAG 
CTGACTAAACAGTAAGTCGACCG 
ATCTGATGCTGACCAACGTC 



Table S9: Codon-tRNA assignments 



tRNA Unmodified Anticodon Modified Anticodon (if known) Cognate Codons Reference 



Leul 


GAG 




cue 


[541 


Leu2 


GAG 




cue, CUU 




[54] 


Leu3 


UAG 


cmo^UAG 


CUA, cue, CUU 




[16] 


Leu4 


CAA 


CmAA 


UUA, UUG 


[54, 61] 


Leu5 


UAA 


cmnm^UmAA 


UUA, UUG 


[54, 61] 


Arg2 


ACG 


ICG 


ecu, CGG, CCA 


54, 61 


Arg3 


CGG 




CGG 




[54] 


Arg4 


UCU 


mnm^UCU 


AGA, ACG 


[<), (il] 


Arg5 


GCU 




ACG 




54] 


Serl 


UGA 


cmo^UGA 


UCA, UCU, UCG 


[54, 61] 


Ser2 


GGA 




UGG 




54] 


Ser3 


GGU 




ACG, AGU 




'54] 


Ser5 


GGA 




UCC, UCU 




'54] 


Prol 


CGG 




CCG 




54] 


Pro2 


GGG 




ccc, ecu 




M] 


Pro3 


UGG 


cmo^UGG 


CCA, ecu, CCG 


[.54, 61] 


Ilel 


GAU 




AUG, AUU 




54] 
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Table S9 (contd.): Codon-tRNA assignments 



tRNA 


Unmodified Anticodon 


Modified Anticodon (if known) 


Cognate Codons 


Reference 


Ile2 


CAU 


k^CAU 


AUA 


[54, 61] 


Glnl 


UUG 


mnm^s^UUG 


CAA, GAG 


[54, 61] 


Gln2 


cue 




GAG 


[54] 


Phe 


GAA 




UUG, UUU 


[51] 
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